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FLOWS DRIVEN BY ROUGH PATHS 



I. BAILLEUL 



("^ ' Abstract. We show in this work how the famihar Taylor formula can be used in a 

Cn I simple way to reprove from scratch the main existence and well-posedness results from 

f— I . rough paths theory, in the context of Holder weak geometric or branched rough paths; 

j^ ' the explosion question, Taylor expansion and Euler estimates are also dealt with. Unlike 
other approaches, we work mainly with flows of maps rather than with paths. We 

Q^ ■ illustrate our approach by proving a well-posedness result for some mean field stochastic 

^NJ ' rough differential equation. 
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1. Introduction 

Since Lyons wrote his groundbreaking article [Ij on rough paths, there has been a con- 
stantly growing industry in understanding the scope of the theory. Besides providing an 
alternative enlighting view on Ito's theory of stochastic integration and some deep asso- 
ciated results, rough paths theory has now invaded the world of numerical simulations, 
stochastic and deterministic partial differential equations and finance, to name but a few 
areas. Despite Lyons' Saint Flour lecture notes [2J, his book [3j with Qian and the im- 
pressive and exhaustive book [3] of Friz and Victoir, rough path still seems to be seen as 
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2 I. BAILLEUL 

a somewhat difficult and teclinical subject wliere algebra and classical analysis meet in an 
intricate way. 

We show in this work how the main existence and uniqueness results of the theory can 
be proved from scratch using only Taylor formula, in the context of dynamics driven by 
Holder weak geometric or branched rough paths. Contrary to Lyons, Friz-Victoir and 
Gubinelli's approach, we work primarily with maps from the state space E to itself rather 
than with E-va\ued paths. Our dynamics on the space of maps cp will have typical form 

(1.1) dip = Vdt + FdX, 

for some driving vector fields V and F = (Vi, . . . , Vi) and a (Holder weak geometric or 
branched p-) rough path X on some time interval [0,T]. Roughly speaking, a rough path 
consists of an E'-valued non-smooth path X together with a number of objects which 
play the role of the missing iterated integrals of X, in the sense that they satisfy the 
same algebraic relations as the iterated integrals of any smooth path, as well as some 
natural size requirements. Given such a rough path and some vector fields V,Vi, . . . ,Vi, 
a family (/ijs)os£s^ts£T of diffeomorphisms of E is constructed from an ordinary differential 
equation. These maps do not form a flow, in the sense that fits is not equal to fifu ° fJ'us for 
allO^s^u^t^T. However, some conditions on the driving vector fields and the rough 
path ensure the existence of a unique flow of maps (</?ts)o^s^t^T close to {fits)o^s^t^T- This 
is what we call a/the solution to the rough differential equation (jl.ip . This approach mixes 
Davie's original approach [5] and Feyel-de la Pradelle's picture [6], [7]. 

All other approaches to rough paths consider paths associated with a point motion as 
the fundamental basic object. Lyons, and later Gubinelli, interpret (jl.ip as an integral 
equation, which requires a suitable notion of rough path integral as a mechanism to attach 
to a rough path Y and some sufficiently regular one form g another rough path L g(y) dY. 
Solving equation (jl.ip then amounts to find a fixed point to an integral equation of the 
form Y = J„ g(Y) dY, where Y is some extension of the original rough path X. See for 
instance [T], [2j or [3], and Friz and Hairer's forthcoming review of Gubinelli's approach. 
Friz and Victoir's approach, developed after Davie's ideas [5], consider M -valued dynamics 
driven by a rough path as a limit in some appropriate topology of elementary dynamics 
generated by controlled ordinary differential equations in which the control converges in 
a rough path sense to some limit rough path. With such a view, no notion of integral is 
needed to define a dynamics; neither is it the case in the approach developped in this work. 

Roughness of a rough path is quantified by a real number p > 2. In order to bring 
the reader progressivey to the most general results, we deal in section [2] with Holder weak 
geometric p-rough paths, with 2 < p < 3, for which only one substitute of iterated integral 
is needed. The case of a general Holder weak geometric p-rough paths, for any p > 2, 
is treated in section [3l while fiows driven by branched rough paths are treated in section 
[5l It comes as a nice feature of our approach that convergence estimates for Euler type 
schemes and Taylor expansion for solution flows come almost for free; this is explained 
in section HI Although designed for bounded vector fields, all the machinery works for 
vector fields with at most linear growth and proves non-explosion of the dynamics in that 
case. The results on fiows are applied in section |6] to give a simple proof of the main 
existence/uniqueness results for the classical rough differential equations on paths driven 
by Holder weak geometric or branched p- rough paths. They are illustrated differently in 
section [7] where we prove a well-posedness result for some large class of mean field stochastic 
rough differential equation. 
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Although everything below works word by word in any Banach space, we have chosen to 
work in a finite dimensional setting where the definition of a 7-Lipschitz vector field, with 
7 non-integer, is well-known, and for which there is no need to pay attention to nice tensor 
norms on tensor powers of the Banach space. The reader acquainted with these notions 
will have no difficulty in putting herself/himself in that setting while the others will be on 
a firm ground. 

A few notations will be used thoughout the text, which we gather here. 

• Fix T > and a non-integer 7 > 1, with integer part [7] and fractional part 
{7}. We say that a vector field W on M'^ is 7-Lipschitz if it is C^ with a {7}- 
Holder continuous differential of order [7], and all its derivative are bounded. Its 
7-Lipschitz norm ||VK||^ is defined as 

[7] 
(1-2) \\W\\, = Yl ||T^(^)|L + ||^""'^||{,} < 00, 

r=0 

wheieW^^' is the differential of order r of PV^ and VF*-™' L 1 stands for the classical 

II ll{7} 

{7}-H61der norm of Ty(W). 

• Given any two C^ vector fields V, W on M we write VW for the vector field x 1— )• 
{DxW){V{x)). Given some sufficiently regular vector fields Vi, . . . ,Ve on M and 
a tuple / = (ii, . . . ,ir) € [1,^]]^, we identify vector fields with derivation operators 
and write V/ for the differential operator / G C^' 1— )• ^1 (Vi2(. . . Vi^/)). Writing 
\y, W] for the bracket of two vector fields, this defines a first order differential 
operator, that is a vector field. Set 

V[i]=[v,„[V,„...,[V,^,„V,^]]... 

• We identify in the sequel L(R^) = R'^ (g) (M^)* and M^ (g, M* with Me{R), via the 
matrix representation of linear maps in the canonical basis. In these terms, given 
any two vectors x, y of M^, the [jk)-t]i component (x (8> yy^ oi x y is x^y^. 

• We use the convention that a*6j stands for the sum ^^ a*6j. 

2. Flows driven by Holder weak geometric p-rough paths (2 < p < 3) 

We explain in this opening section the meaning of equation (jl.ip and how to solve it in 
the simplest case where rough paths are really needed. Recall that a flow on M'^ is a family 
(¥'is)o^s«<T of maps from M'^ to itself such that Lpts = ^tu ° ^us, for all ^ s ^ n ^ t ^ T, 
and (fss = Id, for all s. As a guide to our approach, consider the fiow of diffeomorphisms 
of M'^ generated by the controlled ordinary differential equation 

N 

(2.1) yu = Y,fi{yuW{u) 

i=l 

for ^ u ^ 1, for some smooth vector fields /i, • • • , /at on M'^, bounded and with bounded 



derivatives, and some smooth R -valued control h. The map tpts associates to any x G M 
the value at time t of the solution to (j2.ip started at time s from x. This flow can be 
constructed directly as follows, without involving any integral equation. 

Given < t ^ T, denote by ^xts the maps x ^ x + J^f^i {h'{t) - h\s))Vi{x) from R'^ 
to iself; they depend continuously on (s, t) in the topology of uniform convergence on R*^. 
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They enjoy the following fundamental property which justifies calling them approximate 
solution flows; it is a simple consequence of Taylor formula. The inequality 

(2.2) WfJ'tuO fJ'us- fJ-tsW^^colt - s\'^ 

holds for some positive constant cq, for allO^s^n^t^l. 

Theorem 1. There exists a unique flow {ipts)o(is^t^i on W^ to which one can associate 
two positive constants c and 5 such that 

(2-3) \\(pts -Aitslloo ^ c|t-sp 

holds whenever t—s^6. Moreover, the (pts are the uniform, limits offing := fis„sn-i o • • • o 

Hsiso, where Si = s + ^{t- s). 

The above flow is the flow associated with the ordinary differential equation (|2.ip since 
the latter satisfies condition (12. 3p . We give a complete proof of theorem [T] as it illustrates 
perfectly in an elementary setting our approach to flows driven by rough paths. 

Proof - 1. Uniform controls. We first show by induction that 

(2.4) \\Tl'^,-fits\\^^L\t-s\' 

holds for some positive constant L, for all integers n ^ 1 and all ^ s ^ t ^ T 
sufhciently close. Write for that purpose u = s„+i_2 and note that 

+ IJA'tu o fJ-us — l^ts\\^- 

Using the induction hypothesis, inequality (12. 2p and the fact that the maps (/^ts (•))*_ <i 
are uniformly Lipschitz continuous, with a Lipschitz constant no greater than (l + 
0^(1)), we get 

W'Pts^^ -/^t^oo ^ W'Ptu- lJ'tu\\^+ {l + 05{l))ciu-sf+Co{t-sf 

^ L{t - uf + (1 + os{l))c{u - s)2 + co{t - sf 



<{i^y-i^-"M^)>--> 



To close the induction, it suffices to choose 6 small enough to have the above term 
{•••}< 1, and to choose L big enough to have {• • • }L + cq ^ L. 

Applying the above reasoning to the (R'^xL(M'^))-valued maps x G M'^ i— )■ (^fits{x),DxfJ-ts) 
associated with the ordinary differential equation satisfies by (yu)o^u^i and (-Dxyu)o^u^ii 
where yo = x, it follows from the analogue of (j2.4p that D^fl^ = GiZo Dx^Si+iSi re- 
mains bounded, so the maps JT^"^ ; n ^ 2, i — s ^ (5} are uniformly Lipschitz continu- 
ous, with a Lipschitz constant no greater than 1 + ci(t — s), for some positive constant 
ci- 

2. The existence and uniqueness proofs of the statement of theorem [1] both rely on 
the following elementary identity 

(2.5) 

m 
fmO- ■ -o/l -gmO- ■ -ogi = 2_] idmO- ■ ■og.rn-i+iofm-i - gm°- ' 'O gm-i+l° gm-i]ofm-i-l°- ' "o/l, 

where the gi and fi are maps from M to itself, and where we use the obvious convention 
concerning the summand for the first and last term of the sum. In particular, if all 
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the maps gm ° ■ ■ ■ ° dk are Lipschitz continuous, with a common Lipschitz constant c', 
then 

m 



(2-6) ||/mO--- o/i -gmO---°gi\\^ ^ c'^ll/i 



a) Existence. Set D^ := |0 ^ s ^ t ^ T; t — s ^ 5} and write D5 for the intersection 
of D(5 with the set of dyadic real numbers. Given s = a2~^° and t = 62""° in D5, 
define for n ^ uq 

(n) _ 

where Sj = s + i2~" and s^tn) = ^- Given n ^ no, write 

Ar{n)-1 



Mt 



(n+1) _ ^^ / \ 

ts — U l^Si+iSi+2-"-i o ^g^_^2-"-lsJ 



i=0 



and use (j2.5p with /j = /^s- , ;^s.+2-"-i °/^Si+2-"-isi ^^^d (^j = ^Sj^^s^ and the fact that ah 
the maps /isjv(„)Sjv(„)_i o • • ■ o Hsj^^„.^_^^^Sf^(„^_, are Lipschitz continuous with a common 
Lipschitz constant c', as established in point 1 above, to get by (|2.6p and (|2.2p 



N{n)-1 

II (n+1) (n)\\ ^ I ST^ II II ^ ntl'-,-n\ 

Wts -MLIL^C 2^ ||/^s,+iSi+2-"-lO/^«,+2-"-ls, -/^..+is,||^ ^0(2 ); 

i=0 

so /i*^"-* converges uniformly on B^ to some continuous function 93. We see that Lp 
satisfies inequality (|2.3p on D^ as a consequence of (j2.4p . since the /x *•"■•* are of the form 
^fc(n.) As c/9 is a uniformly continuous function of (s,t) € B^, by (j2.3p . it has a unique 
continuous extension to D^, still denoted by 93. To see that it is a fiow on D^, notice 
that for dyadic times s ^ u ^ t, we have yjj" = i^\^ o ipus , for n big enough; so 
'^ts = ^tu°Vus for such triples of times in B^, hence for all times since (p is continuous. 
The map 99 is easily extended as a flow to the whole of {0 ^ s ^ t ^ T}. 

b) Uniqueness. Let ^ be any flow satisfying contion (|2.3p . With formulas (|2.5p and 
(|2.6p in mind, rewrite (j2.3p under the form -i/^f^ = /ij^ + Oc((i — ■5)^)7 with obvious 
notations. Then 

'^ts = l/js^nSin^i O • • • O -i/;^^^^, = f^52nS2n_i + Oc(2"p") j O • • • O (^flsiSQ + Oc(2"p") j 

where A„ is of the form of the right hand side of (|2.5p . so is bounded above by 0(2""") 
since all the maps /i^jusj"-! °' ' '°^s2"-f+i'52"-<? ^^'^ Lipschitz continuous with a common 
Lipschitz constant. Sending n to infinity shows that 'i/'ts = '^ts- [> 

With the solution fiow to ordinary (controlled) differential equations constructed as 
above, the reader will easily convince herself/himself that the above proof is sufficiently 
robust to justify the following result, which is at the core of our approach to fiows driven 
by rough paths. 

Theorem 2. Let (/a(s,t; •)) _-, at be time dependent regular enough vector fields on 
M.'^, uniformly continuous on {0 ^ s ^ t ^ T}. Let uts{x,A) = (^fits{x), Mts{x, A)^ be the 
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map from R'^ x L(IR'^) to itself which associates to {x,A) the time 1 value of the ordinary 
differential equation 



yu= ^ fa{s,t; yu) 

a=l..N 



a=l..N 

Z-u 

a=l..N 

started from {x,A). Suppose that there exists a > 1 and a positive constant cq such that 

holds for all ^ s ^ t ^ T . Then there exists a unique flow (v^isjo^s^i^i on M'^ to which 
one can associate two positive constants c and 6 such that 

W^ts -/"tslloo ^ c\t- s\° 

holds whenever t — s ^ 6. 

The approach to rough paths driven dynamics developed in this work is a variation on 
this result and its proof. Holder weak geometric p-rough paths (for 2 < p < 3), which 
control the rough differential equation (ll.ip and play the role of h, are defined in section 
12.11 As M -valued objects, they are not regular enough to provide the crucial property (12. 2p 
and define the solution flow as above. The missing bit of information needed to stabilize 
the situation is a substitute of the non-existing iterated integrals f X} dX^ , which provide 
a description of what happened to X during any time interval (s,t). The second part of 
a weak geometric p-rough path plays precisely that role. An approximate solution flows is 
defined from ordinary differential equation involving the vector fields V,Vi, . . . ,Vi and the 
driving rough path X, and studied in sections 12.21 and 12.31 where we prove an analogue of 
inequality (12. 4p . A solution flow to equation (11. ip is then defined in section 12.41 as a fiow 
uniformly well approximated by the approximate solution fiow; theorem [14] provides the 
main result of this section, under the form of a well-posedness result. The general case 
of a Holder weak geometric p- rough path, for any p > 2, is considered in section [3l while 
branched rough paths and their associated dynamics are the subject of section [5l 

2.1. Holder weak geometric p-rough paths. Given a smooth M'^-valued path {xt)o^t^T, 
it is an elementary, though important, fact that the map X^^ : (s, t) i— t- ixt — Xs,f^ /J dxu <8) 

dxr^ G M'' X (1^*^(8) M^) , takes its values in a particular subset of W'- x (M^(g)R'^) specified by 
algebraic conditions. It is precisely these conditions which appear in the definition below. 

Definition 3. Let 2 < p < 3 be given. A Holder weak geometric p-rough path on 

the interval [0,T] is a function X : {s,t) i— )■ {Xts,^ts) from the set {0 ^ s ^ t ^ T} to 
M X (M (g) M ) satisfying the algebraic requirements 

(i) Xtr = Xus + Xtu, for alio ^ s ^u^t ^T, 
(ii) Xts = X„s + ^tu + Xus '^ Xtu, for alio ^ s ^u ^t ^T, 
(iii) the symmetric part o/Xj^ is 2-Xts 'S> Xts, for any ^ s ^ t ^T, 
as well as the analytic conditions 

(a) X is --Holder continuous: \\X\\i := supo<s<j<T - — ^ < oo, 

^ p ^ ^ \t-s\P 

(b) ||X||2 := sup ^^ <oo. 

P 0^s<ts:T \t~s\p 
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The algebraic conditions satisfied by a Holder weak geometric p- rough path justify think- 
ing of its (M.'^ (8) M'^)-part as a kind of iterated integral of X against itself. Condition (ii) is 
nothing but the analogue of the elementary property J J = J J + J J +//> satis- 
fied by any reasonable notion of integral on M such that J = J + / • Note that condition 

(i) on X implies that Xss = for all s. We write X^^ and X^^ for the components in the 
canonical basis of Xts and X^^ respectively. Note that the space of p-rough paths is not 
a vector space as the algebraic conditions on X are not linear. We define a norm on the 
space of Holder weak geometric p-rough paths setting 

II "V II ll^v^ll \/ll\rll 
||A.|| = ||A II 1 V ||A|| 2 , 

P V 

with associated distance function (i(X, Y) = ||X — Y||, although X — Y is not a Holder 
weak geometric p-rough path. The map X^^ associated with a smooth path is a Holder 
weak geometric p- rough path. The term "weak" was coined so as to emphasize the fact 
that there are Holder geometric p-rough paths, as defined above, which are not limit of 
some X^,. in the metric topology associated with d - see for instance [8] or [3]. 

Holder weak geometric p- rough paths appear naturally in a probabilistic context. Let 
B he a. Brownian motion; the random process B = (Bt — Bg, f B^ (^ odBr) is 



* ' O^s^t^T 



almost-surely a Holder weak geometric p-rough path, for any 2 < p < 3. It is called the 
Brownian rough path. It is elementary to prove that E[||B||'^1 is finite for all r ^ 1, using 
the scaling properties of Brownian motion - see for instance [4], chap. 13. More generally, 
Brownian motion can be replaced by any martingale enjoying the same regularity properties 

as Brownian motion. Note that using an Ito integral, ^ = { Bt — Bs, f B^^ dB^ I 

defines a function satisfying all conditions of definition |3] except (iii) since Bf^B^^ is not 
equal to f BrdB^ + J B^dBl for j = k\ it is not a Holder weak geometric p-rough path. 

2.2. Approximate rough difTerential equation. Given 2 < p < 7 ^ 3, let 1/ be a 
globally Lipschitz vector field on R and Vi, . . . ,Vi be 7-Lipschitz vector fields on M . Let 
X be a Holder weak geometric p-rough path. For ^ s ^ t ^ T, denote by ^ts the map 
from M to itself which associates to any x € M the value at time 1 of the well-defined 
and unique solution of the ordinary differential equation 

(2.7) ^yu = {t-s)V{yu) + XiMyu) + \{^ts-\xts®XtsY[Vj,Vk]{yu), ^ u ^ 1, 

with yo = X- This equation is to be understood as the classical ordinary differential 
equation version of equation (II. ip ; the definition of a solution flow to (II. ip given in section 
12.41 will make that point clear. The following property of iits is the main reason for its 
introduction. Write A for the Lipschitz constant of V . 

Proposition 4. Suppose V is bounded and let M be a common upper bound for the 
(7 — 1)-Lipschitz norm ofVi,...,V£ and the supremum norm of V . There exists a positive 
constant c, depending only on M, A, T such that we have for all f £ C^ 

(2.8) fof,,,-\f + {t-s){Vf)+XiiVJ)+^is{VjVkf)} ^c{i + \\W)\\f\h\t-s\h 

Proof - To shorten notations, write a for (t — s), ¥ for Xl^ and c-'^ for ^{X^i — ^Xst ^ 
Xst} • Let ci be a positive constant such that |6| ^ ci||X||jt — s|p and |c| ^ci||X|||t — 
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2 

s\p . Note that there is a constant C2 depending only on M and T such that we have 
for all ^ n ^ 1 

(2.9) |?/„-x| ^ M(c2 + ||X||)|t-s|p. 

Setting 

W := fiiitsix)) - {x + aiVf){x) +b\V,f){x) + ci\[V„V,]f)ix)], 

we have 

W = tU{Vf){Vu) - {Vf){x)]+b^{{Vd){yu) - {V.f){x)} 



+ ci^{{[V,,Vu]f){yu)-{[V,,Vu]f){^)])du, 



(2.10) V f r{aiVVif)iyr) + lf'{Vi>Vif){yr) + ci'{[Vj,Vkm){yr)}drdu 



which we rewrite under the form 

(*)=a/ {{Vf){yu)-{Vf){x)}du + 
Jo 

-1 ru 

lo Jo ' 

c"" I {[[Vj,Vk]f){yu)-{[Vj,Vt,]f){x)]du. 
Jo 

By the Lipschitz continuity of Vf and (|2.9p . there is a positive constant C3 depending 

only on M, A and T such that the first integral above has absolute value no greater 

than C3M(2 + ||X||) ||/||ci|t - s\^^p. As [Vj,Vk] is (7 - 1)-Lipschitz, hence (7 - 2)- 
Lipschitz, the function [Vj, V^]/ is also (7 — 2)-Lipschitz, so the third integral above 

has absolute value no greater than 2M'^ (^cs + ||X||) ||X||||/||ci |t — s|p. The second 
double integral in (j2.10p is of the form 



1 fU 



b'(aAi + b'' / {Vi>Vif){yr)drdu + ci''a,k), 

with \Ai\ ^ 2M'^\\f\\c2, \Cijk\ ^ 2M^\\f\\c2. Also, by the (7 - 2)-Lipschitz character of 
the vector field Vi'Vi, we have by (12. Dh 

\{Vi>Vif){yr) - {V^,V^f){x)\ ^ M^f\\c. (m(c3 + ||X||)|t - s|p)^~'. 

It follows that 

fif^tsix)) - (fix) + a{Vf){x) + b\VJ){x) + ci''{[Vj,Vk]f){x) + -— {V,Vkf){x) 



is no greater than (c+ ||X||''') ||/||c^ \t — s\p for a postive constant c depending only on 
M, A and T. It remains to use the fact that the symmetric part of X^s is 2^ts '^ ^ts 
to see that c^^ [Vj ,Vk] + ^VjVk = x{^ Vj Vk ■ > 

Remark 5. We only use in the above proof the {'j—2)-Lipschitz character of the vector fields 
VjVk rather than their stronger (7 — 1)-Lispchitz character. There is no longer uniqueness 
for the solution of the ordinary differential equation (j2.7p if the Vi are only {'j — 1) -Lipschitz, 
or if the vector field V is only supposed to be p-Lipschitz, for some positive p. However, the 
above proof shows that ()2.8p holds in that case with the time 1 function associated with any 
solution of equation (12. 7p in the role of iHs-, with the exponent ^ replaced &2/min|l + ^, ^|. 
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If the vector field V is C^ and Vi, . . . ,V£ are C^, the classical results on the dependence of 
solutions to ordinary differential equation with respect to some parameters. For any reals 
a, 6*, c^^, the map exp(ay + ¥Vi + c^^\yj, 14]) (x) associated with the differential equation 

^yu = aV{yu) + VViivu) + o''^ [Vj, Vu]{yu). ^ n ^ 1, 
au 

is a continuous function of (a, 6, c) and a continuously differentiable function of (6, c). The 
following basic fact comes as a consequence of the analytic properties of any Holder weak 
geometric p-rough path and the above topology on the set of Holder weak geometric p-rough 
paths. We write ^ instead of ^ts in the proposition below to emphasize its dependence 
on X. 

Proposition 6. • The diffeomorphisms (/^is)n< <t<T ^^^pend continuously on ((s,t),X), 

in the sense that 

max|||/i^-/iy„||oo,||(/i^)"^ - (^™)"^|L} ^ f{{s,t),{u,v)] X,Y) 

for a continuous function f with /((s,t), {s,t) ; X,X) = 0, for all ^ s ^ t ^ T , 
and any Holder weak geometric p-rough path X. 
• ^4/50, all the nts o^t^ uniformly Lipschitz continuous, with Lipschitz constant no 
greater than l-\-c\t — s\^'P , for a constant c = c' (l + HXlp), with a' depending only 
on M and A. 

2.3. Fundamental properties of the approximate flow. The properties of (/its)o^s^t^T 
described in propositions [71 and 1101 below will be crucial to our analysis of rough differential 
equations made in section 12.41 They correspond to the inequalities (|2.2p and (j2.4p used in 
a crucial way in the proof of theorem [TJ As usual, \\h\\oo stands for the supremum norm of 
a map h from W^ to itself. 

Proposition 7. There exists a positive constant c depending only on M, A and T , such 
that 

holds for all0^s<u<t^T. 

Proof - Fix x € R . To shorten notations, write here a for (t — ti), 6* for X^„, c^ for 
X^^, and primed notations for the same quatities with u and t replaced by s and u 
respectively. In those terms, 

IJ'tuilJ'Usix)) = Husix) +aV{fIus{x)) +6Vj(/ins(x)) + c''' {VjVk){nus{x)) +etu{x), 

where ||etM|| ^ c|t — u| p , with the constant c of proposition HI Using (12. 8j) . (12. 9p and 
Taylor formula, we obtain the following elementary estimates, where ci is a constant 
depending only on M, A and T, and for which we only use the (7— 1)-Lipschitz character 
of the vector fields Vi . 



\aV{fius{x)) - aV{x)\ ^ ci(l + ||X||)|n - s^^, 
(2.11) \b%{fius(.x)) - {b'Vi + b¥v,,Vi){x)\ ^ci(l + max(||X||,||Xp)) 
\o>HVjVk){fius{x)) - c^HVjVk){x)\ ^ci(l + ||X||)||X|||tx-s|p. 



3 

u — s\p , 
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We thus have 
(2.12) 

f^tu{l^us{x)) =x + a'V{x) + 6'V,(x) + c'^\VjVk)ix) + aV{x) + b'Vi{x) + b''b\Vi,Vi){x) 

+ c^''{VjVk){x) + etu{x) + 4,(x), 

where 

|4,(x)| ^ 4ci|l + max (||X||, ||Xp)||n - s|? ^ 4ci(l + |iXp)|t - s\p 

and etu(x) ^ C\t — u|p ^ C\t — s\p , for some positive constant C. Noting that we 
have 

a' + a = {u — s) + {t — u) = t — s, 

b + 6* = X^g + Xl^ = Xl^, 

a'^^ + ai^ + yib^ = xil + X^^ + Xi.Xt = Xi^ 

identity (j2.12p and estimate (|2.8p in proposition |4] provide the awaited inequality: 

\fJ-tu{fJ-us{x)) - IJ.ts{x)\ ^ c(l + ||XP)|t - S\P, 

for a positive constant c depending only on M, A and T. [> 

Remark 8. This remark is a sequel to remark\^ as emphasized in the proof, the proof of 
the above result uses only the fact that VjVk is (7 — 2)-Lipschitz, for which it suffices that 
the Vi be ('J — 1)-Lipschitz, so the conclusion of proposition\^ holds in the case where the 
Vi are only supposed to be (7 — 1)-Lipschitz and if the vector field V is only p-Lipschitz, 
with an exponent - in the conclusion replaced by min { 1 + - , - } ; the constant c depends in 
that case of M, p and T. From now on, we shall say that such a constant depends on the 
data of the problem. 

Given a partition vr = {s < si < • • • < s„_i < i} of (s,t), and u G (s,t), denote by 
7r(s, u) and ir{u, t) the partitions of (s, u) and (u, t) induced by vr. If u is one of the Si then 
both cardinals of 7r(s,tt) and TT{u,t) are strictly smaller than n. 

Definition 9. Let e G (0, 1) be given. A partition tt = {s = sq < si < ■ ■ ■ < s„_i < 
Sn = t} of (s,t) is said to be of special type e if we have e ^ s»-^»-i ^ i _ ^^ Jqj. q,// 
i = 1 . . .n — 1. The trivial partition of any interval into the interval itself is also said to 
be of special type e. 

A partition of any interval into sub-intervals of equal length if for instance of special 
type g; the dyadic partition {s < ■ ■ ■ < s + k2~"{t — s) < ■ ■ ■ < t} is of special type 2- 
Given a partition tt = {s = sq < si < ■ ■ ■ < s„_i < Sn = t} of (s, t) of special type e 
and u G {si, . . . , Sn~i}, the induced partitions of the intervals (s, u) and (u, t) are also of 
special type e. Given any partition 7r(s, t) = {sq = s < ■ ■ ■ < Sn-i <t = s.„} of the interval 
(s, t), set 

Proposition 10. Let e G (0, 1) be given. There exists some positive constants cq depending 
on € and the data of the problem, and 6 depending on ||X|| as well, such that for any 
^ s < t ^ T with t — s ^ S, and any partition Tr{s, t) of (s, t) of special type e, we have 

(2.13) ||7i.(,,t) - ^lts\\^ ^ co(i + ||xp)|t - s|p . 
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Proof - Write L for Co(l + HXp), with a constant cq to be determined later. We proceed 
by induction on the number n of sub-intervals of the partition. The case n = 2 is a 
consequence of proposition [71 Suppose the statement has been proved for n ^ 2. Fix 
^ s < t ^ T with t — s ^ 5, and let 7r(s, t) = {sq = s < si < ■ ■ ■ < Sn < Sn+i = t} 
be a partition of (s, t) of special type e, splitting the interval (s, t) into (n + 1) sub- 
intervals. Set u = Spj^, so the two partitions it{s, u) and 7r(u, t) are both of special type 
e. Then 

(2_14) 

^ ||7^7r(M,t) ~/^t«||oo + 11^*" °V'tt{s,u) - fJ-tu o fJ-us\\^ + \\fJ-tu o fJ-us - Mtslloo 

^ L\t-u\i + ('l-Fc(l + ||X|P)(5^)l|u-s|p +c(1 + |lxp)|t-s|^ 

since 7r(u, f ) and it{s, u) are partitions of special type e, and by the induction hypothesis 
together with propositions [6] and [71 Set u — s = a{t — s), with e ^ a ^ 1 — e. The 
above inequality rewrites with shorthand notations 

||7^7r(s,t)-^ts|L ^ {((1-")^ + {l + os{l))ap^L + a^\t-s\i. 
In order to close the induction, we need to choose L and 6 satisfying the condition 
(2.15) a + {(1 - a)p + (1 + os{l))a^L ^ L. 

As e ^ a ^ 1 — e, choose first S small enough to have (1 — a)p + (l + 05(1)) a p < 1, 
then choose L big enough to have a < (l — (1 — a)p — (l -|- os{l))ap)L. The choice 
of (5^/P = 2c7i+||X|h) ' l^^-"^^ foi' instance to L = ^_Q^^x-^/p (l + ||X||'''). [> 

Remark 11. The conclusion of proposition [TU holds under the assumptions that V is p- 

Lipschitz and the Vi are {'^ — 1)-Lipschitz, with an exponent min {l + fj ?} in (|2.13p instead 

of^. 
J p 

It is important to notice here that the statements of propositions l4l [71 and [TOl hold under 
a relaxed hypothesis allowing linear growth of the vector fields V and Vi, under the form 
given below in (|2.16p and (|2.17p . We still suppose that V is Lipschitz continuous and 
the differentials of the Vi of order at most ([p] — l) are bounded and their differentials of 
order ([p] — l) are (7 — [p])-Holder continuous. We talk in that case of (7 — 1)-Lipschitz 
(resp. Lipschitz continuous) vector fields with linear growth. In that setting, the proof of 
proposition [H gives the existence of a positive constant c such that 
(2.16) 

/(/zt,(x))-{/(x)+(t-5)(y/)(x)+X,*,(y,/)(x)+Xi^F,I4/)(x)}|^c{lV|x|}||/||^|t-s|p 

holds for all x € M"^ and any f £ C^ . Similarly, the proof of proposition [TOl gives the 
existence of two positive constants 6 and L such that 

(2.17) \'Pn{s,t)i'') - ^<«(^)l ^ ^{1 V \x\] \t - s\i, 

holds for all x G M"^, all {s,t) G Ds and any partition 7r{s,t) of {s,t) of special type e. 
Inequalities (12.16^ and (I2.17P also hold if we only suppose that V is p-Lipschitz, with the 
exponent ^ replaced by min { 1 + f ) p } ■ 
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Suppose the vector field V is (1 + p)-Lipschitz and the vector fields Vi are j-Lipschitz. 
Recall the constant M appearing implicitly in (j2.13p via cq is an upper bound for the 
(7 — 1)-Lipschitz norm of the 7-Lipschitz vector fields Vi. This gain in the norms is 
useful in that it enables us to apply (j2.13p to the dynamics asscociated with the ordinary 
differential equation satisfied by the differential with respect to x of the solution {yu)o^u^i 
to equation (|2.7p . involving the p-Lipschitz, resp. (7 — 1)-Lipschitz, vector fields {x,A) 1— )• 
{V{x), {D.^V)A), and {x,A) ^ {Vi{x), {a^V{}A) on W^ x L(M'^). Set 

Ats{x) = DxVi = D^fitsi-). 

Then, by remarks [5l [HI and [TTl estimates (j2.16p and (|2.17p hold for Ats] as Agg = Id for 
all ^ s ^ T, it follows in particular from ()2.17p that A^u^f\ = Djl^u^^-^ is bounded, 
independently of the special partition 7r(s, t) and ^ s ^ t ^ T with t — s ^ 8. 

Corollary 12. Fix e > 0. There exists some constants c and 5 depending only on the 

data of the problem and |jX|| such that all the maps {//j^fs t)(') > '""(■S) t) of special type e, ^ 

s ^ t ^ T, t — s ^ 5^ , are Lipschitz continuous, with Lipschitz constant no greater than 
1 

1 + c\t — s\p . 

By proposition [TOl the constant c of the above corollary can be chosen to depend on X 
as a multiple of 1 + ||X|1'''. 

2.4. Flows driven by Holder weak geometric p-rough paths. 

Definition 13. Let 2<p<j<3be given, V be a globally Lipschitz vector field on W^ 
and Vi, . . . ,Vi be C'^ vector fields on M . Let X &e a Holder weak geometric p-rough path. 
Write F for (Vi, . . . ,Ve). A flow ((/3st)os;ss£iCT is said to solve the rough diff"erential 
equation 

(2.18) dip = Vdt + FdX 

if there exists a constant a > 1, independent of^, and two possibly ^-dependent positive 
constants 5 and c, such that 

(2.19) \\ipts-fits\\^^c\t-s\^ 
holds for all ^ s ^t ^T with t — s ^ 6. 

Note that any solution fiow depends continuously on {s,t) in the topology of uniform 
convergence on M by proposition [6] and (I2.19P . 

Theorem 14 (Well-posedness). Given < p < 1 and 2 < p < j !^ 3, let V be a 
{1 + p)-Lipschitz vector field on M. andVi, . . . ,V£ be j-Lipschitz vector fields on M.. Let^ 
be a Holder weak geometric p-rough path. Then the rough differential equation (I2.18P has 
a unique solution flow. It takes its values in the space of uniformly Lipschitz continuous 
homeomorphisms o/M and depends continuously on X. 

Proof of theorem [14] - We follow the proof of theorem [T] step by step, rewriting the 
details for the reader's convenience. Recall 6 > has been fixed in proposition 1101 and 
set Bs := {{s,t) ■,0<sf^t^T, t - s i^ 6} . Write a for min (l + ^, ^). 

a) Existence. Write D^ for the intersection of D5 with the set of dyadic real numbers. 
Given s = 02""" and t = 62""^ j^ B^, define for n ^ no 

(") _ 
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where Sj = s + i2 " and Sjv(n) = i- Given n ^ no, write 

, _,, N(n)-l 

(n+1) _ ^A f ^ 

fJ-ts — U \f^Si+iSi+2-"-^ ° ^J'Si+2-"-^sJ 

i=0 

and use (12. Sp with fi = f^Si+iSi+2-"--^ °/^Si+2-"-isi ^^^ 9i — f^Si+iSi and the fact that ah 
the maps /Usjv(„)Sjv(„)_i o • • • o^sjv(„)->+iSjv(„)-, are Lipschitz continuous with a common 
Lipschitz constant c', by corollarvll2| to get by (j2.6p and proposition [7] 

N{n)-1 
II (n+1) (n)|| , / V^ II II , r,(a-l)n 

1=0 

so ^^"'' converges uniformly on B^ to some continuous function tp. We see that tp 
satisfies inequahty (|2.19p on D^ as a consequence of (j2.13p . As 93 is a uniformly 
continuous function of (s,t) G ID5, by (|2.19p . it has a unique continuous extension to 
D5, still denoted by cp. We see that it is a flow on D^ noticing that for dyadic times 

s ^ n ^ t, we have ipf^ = ipl^ o ipus\ for n big enough; so ipts = ^tu ° 'fus for such 
triples of times in D5, hence for all such triples since ip is continuous. 
We extend (p> to the whole of |0 ^ s ^ t ^ T} by partitioning [0,T] into {O < 
5 < ■ ■ ■ < {k - 1)5 < T i^ k6] and setting, for ^ < m and s G [l6, {i + 1)5) and 
t e [m6,{'m + 1)5), 

fts = 'Ptm5° fmS(m-l)5 o • • • o (p(i+2)S {i+l)S ° V{l+l)5s- 

It is elementary to check that this formula defines a flow; it satisfies (|2.19p with c = 
Co(l + ||Xp), by proposition [TOl 

b) Uniqueness. Let tp be any solution flow to equation (j2.18p . Rewrite (j2.19p under 
the form Tpts = fJ-ts + Oc{\t — s|"), with obvious notations; then we have 



i'ts = V's2nS2n_i O • • • O ^g^^^, = (^^^^^^^n.i + Oc(2 ''" ) j O • • • O {^j^siSo + Oc{2 ' 

(2.20) 

= A'sansan-i o • • • o /Us-^^o + A„ = ^^^ + A„, 

where A„ is of the form of the right hand side of (j2.5p . so is bounded above by 
c'c2~^°'~^^"', since all the maps /^sansan-i o ' ' ' ° /^S2n_^, isgn.^ are Lipschitz continuous 
with a common Lipschitz constant c', by corollary [121 Sending n to infinity shows that 

■0ts = fts- 

c) Properties of the solution flow ip. Arrived at that point, we have proved that 

the rough differential equation (I2.18P has a unique solution, with values in C(M , M ). 

As all the bounds used above are uniform for X in a bounded set of the space of 

Holder weak geometric p-rough paths, and each //("■■* is a continuous function of X, by 

proposition [6l the flow ip depends continuously on ((s,t),X). 

To prove that it actually takes values in the group of homeomorphisms of M , we shall 

see in section [3?T] that (99^ )o^s^t^T can be identified with the solution flow of a rough 

^ ^ ^ <— 

differential equation driven by a rough path X constructed from X. 

As all the n\g are uniformly Lipschitz continuous maps of M , by corollary [121 the 
limit maps cpts have that property as well; the same holds for their inverses. > 
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Remarks 15. (1) Note that one can actually prove that the maps ipts areC^ -diSeoin.oicph.isin.s, 

although the above appproach does not make it elementary - see for instance [?]. 
This is however a direct consequence of the above results if we suppose that V is 
(2 + p)-Lipschitz and the vector fields Vi are (7 + 1)-Lipschitz. 

(2) Convergence speed of Euler schemes. Given any partition it = [s = sq < 
si < ■ ■ ■ < Sn-i < Sn = t\ of {s,t), with mcsh |7r| = max (sj+i — Sj), applying 

i=0..n—l 

(|2.6p with fi = fJ,Sijf.iSi o-nd gi = '^Si+iSi '^^^ ^^^ fact that the maps ipi,a form a flow 
of uniformly Lipchitz continuous maps, with common upper bound d = c(2+ HXp) 
for their Lipschitz constants, we have by (|2.19p 

n— 1 n—1 

«=0 i=0 

^c(2 + ||Xp)'r|7r(s,t)P"\ 

This estimate refines and generalizes the corresponding estimate of Lejay |11| in 
our setting with a drift. See also |12| for a similar result, proved using the heavy 
machinery of sub-Riemannian geometry. 



(3) More general drift. The well-posedness result stated in theorem I4 admits an 
easy generalization to the dynamics 

(2.21) dip = Vdh + fdX, 

where h is an a-Holder continuous real-valued function, with a € (2) l] • ^^^ map 
fits is replaced in that case by the time 1 map p^^ associated with the ordinary 
differential equation 

iju = {h{t) - h{s))V{yu) + X4F,(yJ + ^{Xts - \xts ® ^i.}'' [^„ 14](y«). 
Proceeding as in the proof of proposition [^ one shows that 



f ° i^ts - {/ + (Ht) - Hs)) {vf) + xUvj) + nsiVjVkf)} 



<. 



Cf\t 



provided a > |. Defining a solution flow to the rough differential equation (I2.2ip 
as a flow if for which there exists some positive constants 6 and c such that \\ipts — 
Mtslloo ^ c |t — s\^, for all ^ s ^ t ^ T with t — s ^ 6, we can see as above that 
this equation is well-posed under the hypotheses of theorem\14\ for a > |. 

(4) Link with Stratonovich stochastic differential equations. In a probabilistic 



context where ^ is a realization of the Brownian rough path B = (S,]B), theorem I4 



provides a pathwise construction of the unique strong solution of the Stratonovich 
differential equation 



(2.22) dxt = V{xt)dt + Vi{xt) odBf, 

with xq fixed; identity ()2.19p can be seen as a quantitative version of Ito 's formula. 
Indeed, setting xt = ipto{xo), we have for any C^ function g with compact support 



gixt+,)-gixt)-e{Vg)ixt)-{Bt+,-Bty{Vig)ixt)-nl1^,{VjVkg)ixt) 



^c(l + ||B|me^ 
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by (j2.19p and proposition [^ It follows that 

fi 
'g{xt+e) -g{xt) 



E 



a{xs ; s ^ t) 



Vg + ^Y.^;'g){xt), 



e\0 \ ■" 2 . ^ 

1=1 



since E[|[B|p] is finite. The law of (xt)o^t^T *s thus a solution of the martingale 
problem associated with the above Stratonovich differential equation; the general 
theory of stochastic differential equations ensures that {xt)o^t^T is for any xq the 
unique strong solution of (j2.22p started from xq. 

(5) Linear growth of the vector fields V,Vi. It is not difficult to give a version 



of theorem 14 under the relaxed assumptions of linear growth of the vector fields V 
and Vi given above at the end of section \2.'J[ 

Definition 16. A flow ip is said to solve the rough differential equation dip = 
Vdt + -FdX if there exists a constant a > 1 independent of X and two possibly 
^-dependent positive constants 5 and c such that 

\<Pts{x) - ^xts{x)\ ^ c max{l, |a;|} \t - s\°' 

for all x € M , and all ^ s ^ t ^ T with t — s ^ 6. 

The well-posedness result stated in theorem^I^holds in that case. Its proof follows 
the proof of theorem [I^ by replacing all uniform estimates by pointwise estimates; 
the fact that corollary O remains unchanged is crucial for that purpose. The so- 
lution flow is in that case a continuous function of X for the topology of uniform 
convergence on compact subsets ofM.. 

It follows in particular from that result that no explosion can occur, providing an 
alternative proof of the result obtained by Lejay for p-rough paths, with 2 < p < 3, 
in |llj . 5*66 section\^for the general case. 

(6) Controlled weak geometric p-rough path. It is straightforward to adapt the 
above arguments to the case of a driftless dynamics driven by a controlled weak 
geometric p-rough path. The Holder regularity of X and X are replaced in that 

I V I \^ I 

setting by the requirement that sup — — ^^ < oo and sup ' ^"'^ < oo, where 

the control function ui is positive outisde the diagonal and satisfies the inequality 

uj{r, s) + uj{s, t) ^ uj{r, t) 

for allO^r^s^t^T. An example is uj{s,t) = (f)[t) — (t){s), for an increasing 
continuous function (j), or any control uj such that j,u)1\(a is bounded and bounded 
away from on the time interval [0,T]. The time increments t — s ort — r appearing 
in propositions ^ Q and [121 and in the definition of a special partition need to be 
replaced by oj{s,t) and Lo{r,t) respectively. 

There is little loss of generality in working with controls for which the partitions 
of any interval into sub-intervals of equal length is of special type e, for some positive 
€. Adapting accordingly the definition of a solution flow to the rough differential 
equation dip = fd^, theorem [77] holds in that setting. With the above example of 
control, theorem\T^ (and its generalization theorem \24\) extends the results of Lejay 
[llj . Note that by remark (3), and for oj{s,t) = (j){t) — (p{s) , we can add a drift Vdh 
to the dynamics provided the function h is such that \h{t) — h{s)\ ^ c|(/>(t) — (/>(,s)| , 
with ^ < a ^ 1. 
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The remarks of this paragraph also hold for the case of a flow driven by a general 
Holder weak geometric p-rough paths, as described below. 



3. Flows driven by Holder weak geometric p-rough paths, general case 

As made clear by the definition of a weak geometric p-rough path, for 2 < p < 3, a 
rough path takes values in a very special kind of algebraic structure, which we clarify in 
the following section before defining Holder weak geometric p-rough paths, for any 2 < p. 
We shall then give the same meaning to the equation d<p = Vdt + Fd^, by associating first 
an approximate fiow to V,F = (Vi, . . . , Vi) and X, and by showing that there is a unique 
flow close to the approximate flow. 

3.1. An algebraic interlude: tensor algebra over M^ and free nilpotent Lie group. 

For N € NU{oo}, write T^ for the direct sum ^ (M ) , with the convention that (M ) 

■r=0 

N N rj^\ 

stands for M. Denote by a = © a^ and h = (B by. two generic elements of T) '. The vector 

r=0 r=0 

space T^ is an algebra for the operations 

N 

a + b = Q) {Or + br), 

r=0 

(3.1) 



ab = 4) Cr, with Cr =} flfc (8) br-k 
— n ^ — ' 



fc=0 



It is called the (truncated) tensor algebra of M^ (if N is finite). Looking at the 
coefficients ci and C2, we recover the algebraic conditions (i) and (ii) of definition [3] that a 
weak geometric p-rough path needs to satisfy, for 2 < p < 3. 

The exponential map exp : T^ — )• T^ and the logarithm map log : T^ — )• T^ are 
defined by the usual series 

p.n ( _-! \n 

(3.2) exp(a) = ^-, log(b) = J^ ^^(1 - b)«, 

with the convention a*^ = 1 G M C T^ . Denote by Tr^v : T^ — )• T^ the natural 

projection. We also denote by exp and log the restrictions to T^ of the maps vrjv o exp 

and VTjv ° log respectively. Denote by T^ '' , resp. T^ '' , the elements oq • • • ® cat of 

T| such that ao = 0, resp. ao = 1. Then exp : T^ — )• T^ and log : T^ — )• T^ 
are reciprocal bijections. 

The formula [a, b] = ab — ba defines a Lie bracket on T^ . Define inductively E^ = M , 
considered as a subset of T^°°\ and E'^+^ = [E, E""] C T^°°\ 

Definition 17. • The Lie algebra g^ generated by the E^, . . . ,E^ is called the TV- 

step free nilpotent Lie algebra. 

• As a consequence of Baker-Campbell-Hausdorf-Dynkin formula, the subset ex.p (s^^) 

of T^ is a group for the multiplication operation. It is called the A^-step nilpo- 
tent Lie group on M^ and denoted by G) . 
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Note that the restriction to G^ of the projection map -Kpq : T^ — )• T^ , sending 
ao ® ■ ■ ■ ® Oq to oo © • • • © flp, provides a natural projection TTpg from G^ to G^ for any 

Note also that the algebraic requirements (i)-(iii) satisfied by a weak geometric rough 
path, with 2 < p < 3, are actually equivalent ro requiring that X^s = Xts(B^ts should be an 

(2) 

element of € G) , and that we should have X.^s'^tu = ^ts for all ^ s ^ u ^ t ^ T. As a 

— 1 (2) 

consequence, Xt^ = X^q Xto, so the G\ -valued path (Xto)os£i^T determines the 2-indices 
map X : (n, v) i-^ ^vu- 

This is the key remark for proving that the solution flow (p to the rough differential 
equation ()2.18p is a homeomorphism. Indeed, 

/ (n)\~l -1 -1 

Ws ) = f^s^so ° ■ ■ ■ ° l^s„s„., 

can actually be written as JisnSn~i o • • • o Jlsiso, for the time 1 map Jl associated with 
the rough path X„„ associated with the G\ -valued path Xt_,o. As ^ enjoys the same 
properties as /U, the maps (^Js ) converge uniformly to some continuous map 93^ which 
satisfies by construction Lpts ° ^ts ~ ^'^■ 

3.2. Holder weak geometric p-rough path, general case. The relevance of the alge- 
braic framework provided by the A^-step nilpotent Lie group for the study of smooth paths 
was first noted by Chen in his seminal work |13) . Indeed, for any R^- valued smooth path 
ixs)s^o, the family of iterated integrals 

Xfs := ixt - Xs, / dxs2'S>dxsj^,... , dxsi ® ■ ■ ■ ^ dxsj^ 

\ Js Js J S^Sl^--^Sls[^t 

defines for all ^ s ^ t an element of G^ . It suffices to notice that, as a function of t, 
the function X^ satisfies the differential equation 

|xil = xil0dx„ 

in T^ . Weak geometric rough paths are somehow an abstract version of these objects. 

Definition 18. Let 2 < p. A Holder weak geometric p-rough path on [0, T] is a 

G^'^^^ -valued path X : t € [0, T] h^ 1 © A^^ © A^^ © • • • © a]^^ such that 

I A* I 
(3.3) sup ^^ < oo, 

0^s<t^T \t- s\p 

for all i = 1 . . . [p], where we set X^^ = X^^Xf. We define the norm o/X to be 

\xi 



H.s\ 



i ; 



(3.4) ||X|| = max sup 

i=l...[p] 0^s<t^T k _ s|p 

and a distance d(X, Y) = ||X — Y|[ on the set of Holder weak geometric p-rough path. 

Denote by {ei, . . . ,ee) the canonical basis of R^ C T^ and write for a tuple I 
(ii,... ,v) 

ej/] — Cii ) [642; • • • l^h — 1) ^vJJ • • • ) ana e/ — CiiCi2 ■ ■ ■ Sv' 
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where the above products are in T^ . We also write Y^' for the coordinates of an element 

[pi 
Y of T^ in its canonical basis. 

Given a Holder weak geometric p-rough path X, denote by A = © A^ ® • • • © A^P' its 
logarithm, in Magnus-Chen-Strichartz form |10| . 



/ b] \ [p] 

(3.5) exp K] Y. ^tse^i] = E E ^ts'ei = X,., 

\p] 

for all ^ s ^ t ^ T; it takes values in the finite dimensional Lie algebra g^ . Notice that 
since A is polynomial in X, by formula p.2p . it is a continuous function of X. 

Given a bounded Lipschitz continuous vector field V, and some 7-Lipschitz vector fields 
Vi, . . . ,Vi on M , let pts be the well-defined time 1 map associated with the ordinary 
differential equation 

[p] 

(3.6) yu = {t-s)V{yu) + Y. E K'lVyi^ivu). 

The following analogue of proposition H] is our basic step for studying rough differential 
equations. Recall M is an upper bound for the (7 — 1)-Lipschitz norm of the 7-Lipscliitz 
vector fields V^. 

Proposition 19. Let < p < 1 and 2 < p < 7 < [p] + 1 be given. Let V be a 
(1 + p)-Lipschitz vector field on W , and Vi, . . . ,V£ be j-Lipschitz vector fields on M ; let 
X = 1 © X'^ © ■ ■ ■ © X^P' be a Holder weak geometric p-rough path. Set a = min (l + -,-). 
Then there exists a positive constant c, depending only on M, p, T such that 

[P] 

(3.7) fopt,-[f + {t-s)Vf + Yl E ^tfVif} 

r=l leli/f 

holds for any / € C^. The maps pts depend continuously on ((s, t), X) in uniform topology. 

1 
Proof - The proof rests on the following basic estimate: \yu — x\ ^ c\t — s\p , which holds 

for some positive constant c depending on M,p,T and ||X||. Using first the Lipschitz 

continuity of V and the preceeding inequality, we get 

nl [P] 1-1 

f{pts{x))=f{x) + {t-s) {Vf){yu)du+Yl E ^t'/' {V[i]f){ys,)dsi, 

•^° ri=l hell,£Yl -^^ 

= f(x) + {t-s){Vf){x) + 0,{\t-s\-p) + Y, E ^*^"'W iy[i]f)iys.)dsi. 

ri=iheli,eYi •'^ 

The terms A[^' ^ {y[ii]f){ysi) with ri = [p] are of the form A^^^' ^ (V[/^]/)(x) + Oc{\t — 
s[p) since Vr/. 1/ is in that case (7 — [p] + 1)-Lipschitz, hence (7 — [p])-Lipscliitz. For 



^c(2+iixr) 11/11^ it-sr 



the other terms we can write 



Kf' (%]/)(y.J = AL^-'^ (%]/)(x)+o,(|t-.|?)+AS'^^ Y. E ^^tf" r (%]%]/) (y^.)rf^: 



?-2 = l/2G[l/]'- 
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So we have after that first step 

b] 
f{Mx)) = fix) + (t - s)iVf)ix) + E E ^t'/' (%]/)(^) 

ri=i heli/Yi 

^, 1 T T Jo Jo 

where the above sum is on /i € [[l,£]^'i,/2 S [[1,£]]^^, for ri,r2 bigger than or equal to 
1. Repeating the above analysis (\p] — 1) times gives the result in a straightforward 
way since (I3.5P holds. By construction, the constant c of the statement depends only 
on M, A and T. > 

Remark 20. We only use in the above proof the (7 — 1)-Lipschitz character of the vector 
fields Vj rather than their stronger 'j-Lispchitz character. There is no longer uniqueness for 
the solution of the ordinary differential equation (j3.6p if the Vi are only (7 — 1)-Lipschitz 
or ifV is p-Lipschitz, for some < p < 1. However, the above proof shows that (j3.7p holds 
in that case as well, with the time 1 function associated with any solution of equation (j3.6p 
in the role of pts, and the exponent ^ replaced by min {—^, ^)- 

Once estimate (|3.7p is proved, it is elementary to use the (7 — 1)-Lipschitz character of 
the Vi and Taylor expansion as in the proof of proposition [7] to see that the following holds. 

Proposition 21. There exists a positive constant c depending only on the data of the 
problem, such that 

\\ptu o Pus - Pts\\^ ^ c(2 + ||Xp) \t - s\p 
for alio i^s i^ui^ti^T. 



Using this proposition and the analoguous Lipschitz continuity statement for pts, as 
stated in proposition [6] for fits , it is straightforward to proceed exactly as in the inductive 
proof of proposition [10] to justify the following important statement refining the above 
statement. 

Proposition 22. Let e e (0, 1) be given. 

(1) IfV is p-Lipschitz and the Vi are (7 — 1)-Lipschitz, then there exists some positive 
constants L and S depending only on the d of the problem, ||X|| and e, such that 
for any ^ s < t ^ T with t — s ^ 6, and any partition 7r(s,t) of {s,t) of special 
type e, we have 

(3-8) Wp^^.^t)- pts\\oo^L\t-s\'', 

witha={l + ^,^). 

(2) // V is (1 + p)-Lipschitz and the Vi are ^-Lipschitz, then there exists a constant 

c depending only on the data of the problem, e and ||X|| such that all the maps 

{PTT(s,t){') j '''"(s,^) of special type e, {s,t) ^ D, t — s ^ 6^ , are Lipschitz, with Lips- 

1 
chitz constant no greater than 1 + c\t — s\p . 



A careful track of the constants shows that the constants L and c of the above proposition 
depend on X as a multiple of 1 + |[X||'>', as in proposition [2TJ This remark will be useful 
in section [7] when dealing with some mean field stochastic rough differential equations. 
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Definition 23. Let 2 < p < ■^ < [p] + I be given. Let Vi,... ,Ve be C'^-Lipschitz vector 
fields on M , and X 6e a Holder weak geometric p-rough path. Write F for {Vi, . . . , Vi). 
Let V be a bounded Lipschitz continuous vector field on M . With the above notations, a 
flow {ifts ; ^ s ^ t ^ T) is said to solve the rough differential equation 

(3.9) dip = Vdt + FdX 

if there exists a constant a > 1 independent of X and two possibly ^-dependent positive 
constants 5 and c such that 

(3.10) W^Pts - PtsWoo ^ c\t - s\°- 

holds for all ^ s ^t ^T with t — s ^ 6. 

With proposition [22] at hands, the proof of the well-posedness result theorem [T3] adapts 
word by word and provides the following generalization of theorem [141 

Theorem 24. Suppose V is (1 + p) -Lipschitz and the Vi are j-Lipschitz. Then the rough 
differential equation dip = Vdt + F dX has a unique solution flow; it takes values in the 
space of uniformly Lipschitz continuous homeomorphisms o/M , and depends continuously 
on X. 

The continuity of the solution flow with respect to the driving rough path X has the 
following consequence, which justifies the point of view adopted by Friz and Victoir in 
their works. Suppose the Holder weak geometric p-rough path X is the limit of the 
canonical Holder weak geometric p-rough paths X" associated with a smooth (or Lip- 
schitz continuous) M'^-valued path (a^")os£t^T through the data of its iterated integrals. 
It is elementary to see that the solution flow 99" to the rough differential equation (j3.9p 
with driving rough path X" is the flow associated with the ordinary differential equation 
Vu = y{yu)du + Vi{yu)d{x^y . As ||c^" — ip\\oo = On(l)) from the continuity of the solution 
flow with respect to the driving rough path, the flow ip appears in that case as a limit of 
the elementary flows 93". A Holder weak geometric p-rough path with the above property 
is called a Holder geometric p-rough path; not all Holder weak geometric p-rough path are 
Holder geometric p-rough path [8j. 

The deflnition of p^a makes it clear that the geometric framework associated with the 
rough differential equation (j3.9p is the sub-Riemannian geometry with horizontal spaces 

(n) 

spanned by V,Vi, . . . ,Vi. The identity ipto = lim plr, , gives for instance a construction of 

n+00 

the paths (^iptQ{x)^ _^„ as limits of the horizontal paths {Pto {^))n<:f<:T' ^^^ ^^V ■^ ^ ^'^• 

The proof of theoremll4lalso makes it clear that the constant c appearing in the deflnition 
of a solution flow can be chosen to depend on X as a multiple of 1 -|- HXp. It is straightfor- 
ward to prove a well-posedness result for the rough differential equation d(p = V dh-\- fd^ 
along the above lines of reasoning, for a g-Holder path with min {—^, p) > p + 5- Fo^" ^ 
less regular path h we need to construct some mixted integrals involving both h and X to 
prove such a well-posedness result. Consult [14J or [4j, chap. 9.3, to see how this can be 
done. 

For specialists: note that the construction of the solution flow cpts as a limit of the 
pj" provides a straightforward justiflcation of Friz and Oberhauser's theorem about drift 
induced by perturbed driving signals; theorem 2 in |15| . 
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Remarks 25. (1) Convergence speed of Euler schemes. The proof of theorems 



14 and[2^ shows that for any partition vr(s,t) of special type of the interval {s,t) 



n-l 



(3.11) 



into n sub-intervals, we have 

II II /X"^ii II 

\\Pn{s,t) ~ ftsWoo ^ C / ^ llPsi+iSi ~ ^Si+iSi\\oo 

i=0 

n-l 

i=0 
where m = max (sj+i — Sj) is the mesh of the partition and ci = c(l + HXp) . 

i=0..n—l 



This is essentially Friz-Victoir's estimate |12) in the setting of Holder weak geo- 
metric p-rough paths; their exponent — — is slightly better than a = min (l + -, -) 
though. Compare the above simple proof with the sophisticated proof of [12l| based 
on sub-Riemannian geometric arguments. As noted in remark[T^ (4), it is straight- 
forward to get an analogue result for any dynamics dip = FdX on flows driven by 
a controlled rough path. 

Estimate (j3.1ip provides one of the two ingredients needed to state a universal 
bound on the convergence properties of some empirical estimators for E,^g(^ipto{x)^~\ , 
with X € M"^ fixed and any Lipschitz continuous function g with Lipschitz constant 
no greater than 1. Theorem\2^ below applies in particular to Stratonovich differ- 
ential equations driven by Brownian motion; it improves upon Malrieu and Talay's 
results [16\ . Write ip^ and'fi^,^^ to emphasize the dependence of these maps on 

X; for an independent and identically distributed sequence (X*)jj.i, write 'ff , ^ for 

P-K{s,t)- 

Theorem 26. Suppose the conditions of definition \2^ on V and F hold. Let 
(X*)jj.i be an identically distributed sequence of independent random Holder weak 
geometric p-rough paths on the time interval [0,T]. Suppose 5 := E[exp (a||X|p)] 
is finite for some positive constant a. Then there exists some positive constants 
o, h, ci, C2 depending only on M, A and T such that we have 

i=l 

for any partition vr(0, t) of (0, t) of special type with mesh m and any Lipschitz 
continuous function g, with Lipschitz constant no greater than 1. 

Proof - Use inequality (j3.1ip and the Lipschitz character of g to get a straight- 
forward upper bound for P ( E 



|5(p?(o,t)(^)) -9Wmix))\\ > ij. To deal with 
W := IP(|FEiIi5(p;(o,t)W) -IE[<?(p^(0,t)(x))]| > i), we use the fact that 
p^/Q jN is a locally Lipschitz continuous function of X, with Lipschitz constant 
no greater than a constant times R on the ball of center and radius R, inde- 
pendently of the special partition 7r(0,i) of [0,t]. (This constant depends only 
on M, X and T .) Write Y for p^fQ t~)(^)- ^^ have for all ^ u ^ a 



(•)^e-^''{E[e 



K/(n-E[/(y)]] '^ 
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,^u(fiY)-ElfiY)] 



where F(i>) := E 



e 



is finite on [0, a], because Jensen's inequality 



and the fact that ||X|| has a Gaussian moment imply 



F{v) ^ E 



A9iy)-9{Y')) 



C E 



Qmax(!|X||,||X'||) ||X-X' 



< oo 



for Y' independent ofY, with the same law as Y. As logF is strictly convex, 
logF(i^) ^ mill (ciz^^, ^i^), where ci is the variance of f(Y). It follows that 

{-k) ^ exp ( — i/A^(|— min(cii^, — )) I, from which the result follows by optimizing 

over z^ G [0, a]. [> 

Note that Lemaire-Menozzi [17J and Frikha-Menozzi [18j 's results on the Gauss- 
ian concentration property of the elementary Euler scheme for some ltd stochastic 
differential equation suggest that the exponential concentration property of the Mil- 
stein type estimator -^ X^j=i ^(/'^(o t~)('^)) f^'''^[9{^to{x))] can probably be improved 
if some additional mild information on the driving rough path is available. Theorem 
[2B applies to all dynamics driven by any Holder weak geometric p-rough paths with 
a Gaussian moment though. 

(2) Linear growth of the vector fields V,Vi. Proceeding exactly as in remark [731 
(3), there is no difficulty in relaxing the boundedness hypothesis on the vector fields 
V,Vi in theoremWA adopting definition\16\for the notion of a solution flow. 



(3) Non-explosion. It follows in particular that there is no explosion of the solution 
flow in the above linear growth setting. This generalizes the results of |llj and |19j 
to the case of a dynamics with a drift. Note that these two works deal with paths 
rather than with flows of maps and that |19| works with a more general notion of 
rough path than the notion adopted here. 

(4) Full rough diff'erential equation. As first noticed by A. Davie [5\ in a rough 
paths context, a vector field W on M naturally lifts to a vector field on the [p]-step 
nilpotent Lie group g'J^^^ C 0l=o(J^'^)®*) setting W{a) = a(0©P^(ai)eOe- • -00), 
for a = (1 © oi © • • • ® aui) G Gd ■ The vector fields Vi are unbounded 7- 
Lipschitz with at most linear growth if the Vi are 'j-Lipschitz; as they are tangent 
to Gf , their flows leave GJ^ invariant. Theorem [^ and remark (2) above 
apply, providing a flow on GJ^ . The map (^ipsQ(ld)^^ipto{U))^^ _^_^„ is then a 
Holder weak geometric p-rough path which can be used to construct new dynamics 
by solving some rough differential equations - see pi|, Chap. 10. 4, for a different 
approach to that question and more material on it. 

(5) Infinite dimension. All the above results, and their proofs, hold in an infinite 
dimensional setting where M is replaced by any Banach space. The appropriate 
notion of a-Lipschitz vector field, as described in [T], [2] or [3], needs to be used. 
Good norms on tensor products also simplify matters when working with the full 
rough differential equation in a Banach space setting - see the above references. This 
framework is not sufficient though to deal with rough partial differential equations 
where the drift V of section would be a densely defined (potentially non-linear) 
unbounded vector field. 

(6) Time inhomogeneous dynamics. The above result have a straightforward gener- 
alization for a time dependent drift V{s,x) which is a bounded Lipschitz continuous 
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function of its two arguments and vector fields Vi{s,-) which are ^-Lipschitz with 
respect to the space variable and Lipschitz continuous with respect to time. The map 
pts is defined in that case as the time 1 map associated with the ordinary differential 
equation 

[p] 
yu = {t-s)V{s,yu) + ^ ^ Aj/V[j](s,y„). 

r=l/G[l/r 

In particular, inequality (|3.7p holds in that case, with Vs{x) andVj{s,x) in place of 
V{x) and Vi{x). This framework will be useful in section\^ for the study of some 
simple mean field stochastic rough differential equation. 

(7) Backward dynamics. Setting Xj = "Kx-t defines a Holder weak geometric p- 
rough path; as noticed in section \37T\ the inverse of the solution flow ip to the rough 
differential equation dip = Vdt + FdX is the solution flow to the rough differential 
equation d^ = —Vdt + FdX.. 

(8) Flows driven by (non-weak geometric) Holder p-rough paths. Note that the 
only place where we have used the fact that X takes values in the free nilpotent Lie 
group G[ is in writing its logarithm as a sum of terms involving only the brackets 
of the basis vectors ei,...,e£. The above reasonnings are totally independent of 
that fact. Recall all elements ofT^ are invertible. We define a Holder p -rough 

path as a Tf -valued path X such that the Holder type conditions p.3p holds. 
One can think of Brownian motion, together with its Ito first iterated integral, as 
a typical (non-weak geometric) p-rough path, for 2 < p < 3. Set log Xts = A^g = 

X^r=i 'l2ieUiV ^"ts ^l> '^'^^'^ define in that case pts as the time one map associated 
with the ordinary differential equation 

[p] 

r=l/e[[l/]'- 

The above proofs show that theorem [^ holds in that case as well, generalizing 
Lejay's results in [TIJ, who only deals with the case 2 < p < 3, without a drift. The 
point of working with weak geometric rough paths is the fact that only brackets of 
vector fields appear in the definition of pts- These brackets are intrinsically defined 
vector fields which do not depend on any choice of a coordinate system, so are 
well-defined on a manifold. This is not the case of the vector fields Vj which are 
needed in the above definition of pts and are meaningless in a manifold setting. 
This fact justifies calling G^ -valued rough paths (weak) geometric rough paths. 
See the remark closing section\^for a similar situation, in the setting of branched 
rough paths. 

4. Taylor expansion of solution flows to rough differential equations 

It is a nice feature of our approach that Taylor expansion of solution flows to rough 
differential equations come almost for free. As in the classical case, the smoother the 
vector fields Vi are, the better we can describe the solution flow to the driftless equation 

dip = FdX^. This takes the form of a refined version of the inequality \\ipts — Pts\\oo ^ 

I- 
c\t — s\p . 
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Let 2 < p and an integer N ^ [p] + 1 he given. Lyons devised in his original theory of 
rough paths [1] a fundamental mechanism which enables in particular to extend uniquely 

any Holder weak geometric p- rough path to a very special G^ -valued map and provides 
the missing higher order iterated integrals needed to write a Taylor expansion of the solu- 
tion flow to the equation dip = FdX.. We state it here under the form we need - see [7] for 
a simple proof of a refined version of Lyons' extension theorem. 

Theorem 27 (Lyons' extension theorem). Let 2 < p be given and X 6e a Holder weak geo- 
metric p-rough path. Let N ^ [p] + l be an integer. Then there exists a unique (G\ -valued) 
Holder weak geometric N -rough path path Y extending X in the sense that iT^^j^y^ts = '^ts, 

for all ^ s ^ t ^ T . The map Y is called the lift ofK. to G^ ; it is a continuous function 
ofX. 

Let 2 < p, an integer N ^ [p] + 1 and < a < 1 be given; set 7 = A^ -|- a. Suppose that 
the vector fields Vi are 7-Lipschitz. Given a Holder weak geometric p-rough path X, let Y 
its lift of X to Gf^^ and A = e A^ e A^ e • • • e A^ be the logarithm of its Hft, in the 
sense that exp Aj^ = Y^^, for all ^ s ^ t ^ T. 

Denote by p^ the diffeomorphism of M'^ which associates to any x € M"^ the value at 
time 1 of the well-defined and unique solution of the ordinary differential equation 

, N 

r'=l/G[l,^r 

with 2/0 = X. The proof of proposition 1191 shows that the following awaited estimate holds. 

Proposition 28. Let2 <p and [p]-|-l^iV<iV-|-a<A^+l he given, with N integer. 
LetVi,...,Vi be {N + a)-Lipschitz vector fields onW^. Lei X = 1 X^ • • • A^ be 

a Holder weak geometric p-rough path with lift Y to G^ . Then there exists a positive 
constant c, depending only on M,X,T and |[X|[ and f ^ C^ such that 

N 

(4.1) \\f°pts-{f + Y. E "^ts'vif) 

^=i/eli,€]'- 
holds for all f € C^. The maps p^l. depend continuously on ((s,t),X) in uniform topology. 

As in section [3?2l this fundamental estimate implies a well-posedness result. 
Theorem 29. The following holds under the hypotheses of proposition W^ 

• There exists a unique flow {'Pts)o<,s^t<,T on M for which there are two positive 
constants 5 and c such that 

(4.2) hfs-pfslL^clt-sl"^ 

holds for all ^ s ^ t ^T with t — s ^ 6. 

• The flow takes its values in the space of uniformly Lipschitz continuous homeomor- 
phisms ofW^, and depends continuously on ((s,t),X"). 

Let now 7 e (p, [p] -|- l). Since ||/9^ — pts\\ ^ c\t — s\p , the identity 

lkiI-pts|L ^ c\t-s\p 

holds for all ^ s ^ t ^ T sufficiently close, so if^ is a solution flow to the rough 
differential equation dip = FdX., where the Vi are considered as 7-Lipschitz vector fields. 



N + a 

^ c\t — s\ P 

00 
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As it has a unique solution by theorem [2H it needs to coincide with the solution flow ip 
constructed in the proof of theorem [2H 

Corollary 30 (Euler estimates/Taylor expansion). Let 2 < p, an integer [p] + 1 ^ iV 
and < a < 1 be given. Suppose the vector fields Vi, . . . ,Vi are (N + a)-Lipschitz. Let 
X 6e o Holder weak geometric p-rough path. Denote by ip the unique solution flow to the 
driftless rough differential equation dtp = FdX.. Let Y be the lift ofS^ to G\ . Then there 
exists two positive constants 5 and c such that 

II Afll I | ^+" 

Wt^- PtsWoo ^C\t-S\ P 

holds for all ^ s ^t ^T with t — s ^ 6. We have in particular, 

N 

(4.3) \\pts-{ld + Y^ Y. Y^s'Vi} 

for all ^ s ^ t ^T with t — s ^ 6. 

Remarks 31. (1) In a probabilistic context where X is the random realization of the 

k I 

Brownian rough path, the Y^^ coincide almost-surely with the iterated Stratonovich 
integrals L odB]^^ ® ■ ■ ■^odBl'i, for I = (ii, . . . , i^), and (14. 3p is a pathwise version 
of Azencott's celebrated stochastic Taylor formula - see for instance |20) . |21| and 

El. 



N+a 

^ c\t — s\ P 

oo 



(2) Friz and Victoir |12j proved similar estimates by a clever use of geodesic approxi- 
mation in the free nilpotent Lie group G^ . 

(3) Most results of Baudoin from |23| and |24| on the geometry of hypoelliptic diffusions 
can be recovered from estimate (j4.3p without resorting to Chen- Strichartz' formula, 
as done in these works. 

(4) To deal with rough differential equation with a drift we would need to introduce the 
cross iterated integrals between t and X, on which subject we have nothing new to 
add; see for instance [3j, Chap. 3.3.3, for their definition and construction. It is 
straightforward to adapt the above reasoning to that setting. 

5. Flows driven by branched rough paths 
Let us say that a family {i^tsjQ^ <t<T °^ i^aps is an approximate flow if the inequality 



loo 



holds for some positive constants c and a > 1, and all s ^ n ^ t sufficiently close. It should 
be clear to the reader by now that the essential ingredient in our construction recipe of a 
flow on M is an approximate flow on M x L(]R ) obtained as the time 1 map associated 
to an ordinary differential equation, as emphasized in theorem [2j We have seen above that 
such approximate flows can be constructed from Holder weak geometric p-rough paths. The 
appearance in that context of the iV-step nilpotent Lie group comes from the fact that the 
family of iterated integrals of a smooth path need to satisfy the usual rules of calculus 
|13| . which lead in a stochastic setting to Stratonovich-like dynamics. It is however clear 
that a wealth of situations cannot be modelled by such dynamical schemes which somehow 
anticipate over the future, and that Ito-like dynamics are sometimes more appropriate. It 
is possible to extend the machinery of flows driven by rough paths, as descrived above, 
to this non-geometrical setting by using the notion of branched rough paths introduced by 
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Gubinelli in [25| . This section is devoted to explaining how this can be done in a simple 
way. 

Consider an ordinary differential equation 



controlled by some smooth M -valued path h. Given any tuple / = {ii, . . . ,ik) € [!,£]], 



write H^^ for J^^^ <—<s <t'^^s\ • • • ^^\- '^^^ point of departure towards branched rough 
paths is the rewriting of Taylor formula 

N 



under the form 

(5.1) xt = Xs + Y,HlVr{xs) + 0{\t 



^r+^i 



where all the multiple derivatives Vi{x) = T^^ (V^j ( • • • Vi^_-^Vi^) . . . )(x) have been ex- 
panded entirely so as to write each of them as a sum of terms of the form ^^ D'^Vi^ (• • • )j 
where D'^Vij, stands for the a*^ derivative of Vij,. We have for instance ^1(^2(^3)) (2;) = 
(-Da;V3)((ViV2)(x)) -|- {D'^Vs)(Vi{x),V2{x)Y The best way to represent the resulting sum 
is to index the sum by labelled rooted trees, as first discovered by Cayley, which is what 
the above sum over r is. The coefficients H^g are then linear combinations of the initial 
coefficients H^^ . As the family of iterated integrals of h, the family {H'^)t of these coef- 
ficients happens to take values in a very special kind of algebraic structure: the dual of a 
Hopf algebra. Section 15.11 provides an elementary presentation of this structure. Starting 
directly from such objects, subject to some size requirements similar to those satisfied by 
a Holder weak geometric p-rough path, defines a branched rough path. We shall see in 
section 15.21 that it is easy to construct an approximate flow from a branched path and 
vector fields, and define as above a notion of solution flow to a differential equation driven 
by a branched rough path. 

5.1. Branched rough paths. Let 7" stand for the set of (possibly empty, unordered) 
rooted labelled trees, with labels in {1, . . . ,^}. Given ti and T2 in T and a G {1, ...,£}, 
denote by [rir2]a the element of T obtained by attaching the two trees ti and T2 to a new 
root with label a. Any element of 7" can be constructed in this way starting from the 
empty tree. We denote by |r| the number of vertices of a tree r and deflne its symmetry 
factor (t{t) recursively 

a([<^..r,"'=],)=nl!...nfc!a(rl)"l...a(rfc)"^ 

where ri, . . . ,rfc are distinct trees with respective multiplicities ni, . . . ,nfc. Set Xf^ = 1 
and V0 = 0. and deflne recursively, for a G [1, ^J and ri, . . . , r^ in T, 



(5.2) 



v.. = Va, F[.,...,,]„ = \ {D-Va) (K. , . . . , KJ 

Cr[[Tl...Tn\a) 
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In these terms, (|5.ip rewrites 

Xi — Xs 



rer, |r|s;Af 



Formal series of the form 



reT 

for y G M°' and real valued coefficients X'^, are called B-series, after Butcher [26] who 
first considered them implicitly in his work on order conditions for Runge-Kutta methods. 
Their systematic study started with the work of E. Hairer and Wanner |27) . who established 
a fundamental result about composition of B-series. 

Theorem 32 (E. Hairer- Wanner). Let X and Y be two real valued maps on T, with 
X® = 1. Then B{Y ,B{lL,y)^ is again a B-series, say B(Y -k'K,y), where Y-k'K is 
explicitly computable in terms of^ andY'. 

This map (X, Y) i— > Y * X is the reason for the appearance in this context of the above 
mentionned Hopf algebra structure. We follow the forthcoming enlighting work [28j to give 
a short description of this algebra and define branched rough paths. 

Let (T-L, •) stand for the set of commuting real-valued polynomials with indeterminates 
the elements of 7" and polynomial multiplication operation; equip the tensor product 71®% 
with the induced product (a ® h){c ® d) — (oc) ® (bd). We define a coproduct A on "H as 
follows. Given a labelled rooted tree r, denote by Sub(r) the set of subtrees of r with the 
same root as r. Given such a subtree s, we obtain a collection ri, . . . , r„ of labelled rooted 
trees by removing s and all the adjacent edges to s from r. Write r\s for the monomial 
Ti . . . r„ . One defines A : T —^ H (S) Ti hy the formula 

At= Yl (As)®s, 
seSub{T) 

and extend it to H by linearity and by requiring that A(ri . . . r„) = A(ri) . . . A(rn). This 
coproduct is coassociative ((A Id)A = (Id (8) A)A) as should be in any Hopf algebra. 
Given a real- valued map Y on T, we extend it into an element of the dual space to Ti setting 
yTi...r„ _ YYl^^ Y[ , for a monomial ri . . . r„, and by linearity. The map (X, Y) i— > Y * X 
in theorem [32] is given by 

(Y * X)^ = Y ^^^'^' = (Y X) (At) , 

seSub(r) 

for a labelled rooted tree r. The third ingredient needed to define the Hopf algebra structure 
of % is an antipode (3, that is a map <3 : Ti ^^ Ti inverting A in the sense that M(Id(8)6) A = 
M(6 (8i Id) A = Id, where M stands for the multiplication map M(a (8> 6) = ab € T-L. An 
explicit formula for & was first obtained in [29j; see |30j for a simple and enlighting proof. 

Denote by 7i* the dual space of Ti and by = ri . . . r„, with Tj G T, a generic element 
of the canonical basis of Ti. Write (6*) for the dual canonical basis and use bold letters 
to denote generic elements of T-L*; we use the notation [a, 6) or a^ for the duality relation 
between "H* and H. With these notations, for an element a of Ti* of the special form 
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a = X^^/gT-a^ (''"')*) ^^ have for instance 

a*a = ^(a0a)(Ae)r = ^(a «) a)(AT) r* 



reT 

where the inside sum is over the set of subtrees ri of r such that r\ri is a subtree T2 of r. 
More generally, we have for such an a G 7^* 



z^« •••« 



(5.3) a- = Y, 

reT 

where the inside sum is over the set of disjoint rooted subtrees ri, . . . , r„ of r, with respec- 
tive roots ai, . . . ,an such that Oj+i is a descendant of Oj in r, for all i = l..n — 1, and any 
node of r is in one of the subtrees Tj. 

We define a Lie bracket on Ti* setting [a, b]j, = a * b — b • a, for which the real vector 
space spanned by the dual trees is a Lie algebra, say g. The exponential map exp^^ : Ti* — > 
Ti* and the logarithm map log^ : Ti* —?■ Ti* are defined by the usual series 

(1 - a)*" 



exp^(a) = J^ -T' log*(^) = Y^ ' 



n\ '^ — ' n 

with the convention a**^ = 1, where 1 denote the dual of the empty tree. Formula (j5.3p 
justifies the convergence of the these sums, in the sense that (expj^(a),0) and (log^(a),0) 
are actually finite sums for any monomial 6. 

It is clear from the above formula for a*" that we have for a nonempty tree r 

\r\ 

(5.4) (exp.(a),T)=j;-5]a-...a-, 

n=l 

with the same inside sum as in (|5.3p . Setting T-Lq = |a G Ti* ; (a, 1) = O} and T-Lq = |a G 
Ti* ; (a, 1) = 1}, one can see that exp^ : T-Lq -^ Til and log^ : T-L^ -^ T-Lq are reciprocal 
diffeomorphisms; see for instance [31] . As a consequence, one can prove as in the tensor 
product algebra case, [32], theorem 3.2, the following result. 

Theorem 33. The pair (^exp^{g),-k) is a group. 

Write 'HJi^\ for the vector space spanned by the monomials (ti . . . t„)*, with Y17=i l"^*! ~ 
k, and denote by Trjsf the natural projection from Ti* to the quotient space 7i*\ ©fc>7v+i '^Ik) ■ 
The image of exp^(0) by ttat is then diffeomorphic to the real vector space (TJy) spanned 
by TTj\f(T'*); it is a Lie group when equipped with the operation vr^v ° *, which we still 
denote by -k. Write Qj, for that Lie group; it plays for branched rough paths the role that 
the N-step nilpotent Lie group plays for Holder weak geometric p-rough paths. Denote by 
•* the element of the dual canonical basis dual to •„. It is elementary to check that the 
'H*-valued function Xts = (-^ts)TeT, |t|s£7V associated with a smooth M^-valued path h, as 
defined in (j5.2p . considered as a function of t, satisfies in ■n]\j{T-L*) the ordinary differential 
equation 



dXts = Y,2its*{*ldht), 
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which imphes that Xts is an element of Q^ for all t ^ s. We follow Gubinelli [25j . and 
[28j . in defining a branched rough path as follows. 

Definition 34. Let 2 < p. A branched p-rough path on [0,r] is a Q^'^^'' -valued path 
X such that 

I Yd I 
sup ^ < OO, 

for all monomials 9 = ti . . .tj^ with \9\ = X]j=i \Ti\ ^ \p\, where we set 'K.tg = X~^ -k X-t, 

d X^^ stands for the inverse of 
We define the norm ofJ^ to be 



and X^ ^ stands for the inverse o/X^ in the group {Gg 



/'c: r\ ll"V"ll I is I 

5.5 X = max sup ^r, 

9;\d\<S:[p] 0^s<t^T U_ gl^ 

where 6 runs over the set of elements of the canonical basis of %; we define a distance 
(i(X, Y) = ||X — Y|j on the set of branched p-rough path. 

It is elementary to check that the inverse a~^ of any element a of ^^ is given by 
the formula a~^ = S*a, where &* is the dual map in H* of the antipode map in Ti. The 
framework of branched rough paths extends the framework of Holder weak geometric rough 
paths, in so far as any weak geometric rough path can be identified with a branched rough 
path, as proved by M. Hairer and Kelly in |28| . The converse also holds true in some sense, 
so the two notions appear equivalent for practical purposes. 

5.2. Flows driven by branched rough paths. It is straighforward to use the machinery 
of approximate flows, as summed up in theorem [2l to produce a flow from a branched rough 
path and some vector fields. Let then 2 < p < 7 < [p] + 1 be given, together with some 
7-Lipschitz vector fields Vi,... ,V£ on M'^, a bounded Lipschitz continuous vector field V 
on M'^ and a branched rough path X. Let A stand for its logarithm, with values in (TJ^), 

exp, I ^ ALr*| =Y^xt0* = Xts. 
,t£Tn / e 

Let pts be the well-defined time 1 map associated with the ordinary differential equation 
(5.6) yu = {t- s)V{yu) + Yl ^tsMVu), 

where the vector fields Vr are defined in (|5.2p . The following analogue of proposition [12] 
holds true under the above hypotheses on the vector fields V,Vi, . . . ,Vi. Recall M is an 
upper bound for the (7 — 1)-Lipschitz norm of the 7-Lipschitz vector fields Vi, and A is an 
upper bound for the Lipschitz constant of V. 

Proposition 35. There exists a constant c, depending only on the data of the problem, 
such that 



fopts-{f + {t-s){vf)+ Y, xi^iyrf)} 



^c(2+iixir)ii/ii^it-s| 



holds for any f & C^ . 
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Proof - The proof of this resuh is identical to the proof of proposition \T9\ and starts 
with the identity 

f{p,,{x))=f{x) + {t-s){Vf){x) + 0,{\t-s\'')+ Yl ^'tl I {Vrj){ys,)dsi 

In Kb] 

= f{x) + {t-s){Vf){x) + 0,{\t-s[^)+ Y. A?s{Vrj){x) 

In Kb] 



+ Yl ^SAS / {Vr,Vrj){ys,)dS2dsi 

^r 1 -^0 



Inl+lrzKb] 
Writing r2 = [rj • • • T|*]a, we have 

n 

(5.7) K,K, = (I)"+V,)(Kl,v;l,...,K,")+5](Z?"K)(K-l,•••,KlK■|'•••'^-?)• 
The result follows by induction from that identity and identity (j5.4p satisfied by any 
element of the group Ql . [> 



As in section 13.21 the above result still holds true with any solution of equation (j5.6p 
if the vector fields Vi are only supposed to the (7 — 1)-Lispchitz and the vector field V is 
supposed to be /9-Lipschitz, for some < p < 1. The crucial step in constructing a flow 
from the family of maps pts is the following counterpart of E.Hairer and Wanner theorem 
on B-series. 

Proposition 36. There exists a constant c depending only on the data of the problem 
such that 

\\ptu° Pus -Ptslloo ^ '^(^+ Il^ir)l*-'SP 

holds for allO^s^u^t^T. 

Proof - It is a direct consequence of proposition [35l Indeed, as we have for any t gT 

Vr{pus{x))=Vr{x) + {u-s){VVr){x)+ ^ Xl',{Vr'Vr){x) + Oc {\u - s\v) 

k'l=£b] 
with C = c(2 + ||X||T), it follows that 

PtuiPusix)) =X+{t- S)V{X) + Y ^usVrix) + Y ^tX' (x) 

+ Y Xl^K'siVr'Vr){x) + Oc{\t-r\i). 

|T| + |r'Kb] 

The result follows in a straightforward way from proposition [35] and identity (|5.7p , 
since X^^ • Xt„ = Xf^. > 

Starting from that point, the construction of a flow from the approximate flow {pts)o^s^t^T 
follows line by line the construction detailed in sections [2] and [3] for flows driven by Holder 
weak geometretric p-rough paths. We adopt deflnition [23] for the solution flow of a rough 
differential equation 

(5.8) dip = Vdt + FdX 

driven by a branched rough path X. 
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Theorem 37. Assume that V is {1 + p)-Lipschitz and the Vi are ^-Lipschitz. Then, the 
rough differential equation (j5.8p has a unique solution flow; it takes values in the space of 
uniformly Lipschitz continuous homeomorphisms ofM.'^ and depends continuously on X. 

All the remarks of section [3] hold in that setting, except the remark about the full rough 
differential equation. The analogue of the latter would lead us to talk about controlled 
rough paths |33) . about which we have nothing new to add. The results of section |4] 
on Taylor expansion of solution flows also hold in that setting, as an analogue of Lyons' 
extension theorem holds for branched rough paths [25]. Such expansion results seem to be 
new in that setting. 

Although the framework of branched rough paths extends the framework of weak geo- 
metric rough paths, the construction of a flow driven by a branched rough path is intrinsi- 
cally linked with a choice of coordinates on M , in so far as the definition of the vector fields 
Vr depends on that choice through the derivatives D'^Va, for n ^ 1. So it is not possible 
to define flows driven by a branched rough path on a manifold without the introduction of 
another ingredient. A linear connection on the tangent space can play that role, and there 
is no difficulty in extending the above results to that setting if the manifold is compact. On 
the other hand, the definition of a solution fiow to a rough differential equation driven by 
a weak geometric rough path extends readily to a manifold setting as it only involves the 
basic vector fields Vi and their commutators, which are intrinsically defined vector fields. 

6. Paths driven by rough paths 

The results of the preceding sections provide an easy proof of existence and well- 
posedness for point dynamics driven by a Holder weak geometric or branched p-rough 
path. We adopt the notations of the preceeding sections. 

Definition 38. Let 2 < p < ^ ^ [p] + l be given. Let Vi,. . . ,Vi be C'^ -Lipschitz vector fields 
on W^, and X &e a Holder weak geometric or branched p-rough path on the time interval 
[0,T]. Let V be a bounded Lipschitz continuous vector field. A path (zg ; ^ s ^ T) is 
said to solve the rough differential equation 

(6.1) dz = Vdt + FdX 

if there exists a constant a > 1 independent of X, and two possibly ^-dependent positive 
constants 5 and c, such that 

(6.2) \zt - pts(zs)\ ^ c\t - s^" 
holds for all ^ s ^ t ^ T with t — s ^ 6. 

For (7 — 1)-Lipschitz vector fields Vi and a p-Lipschitz vector field V , the differential 
equations (j3.6p and (|5.6p are not well-posed; in that case, pts{x) stands for the time 1 value 
of an arbitrarily chosen solution started from x. Note that a solution path is necessarily 
continuous by (j6.2p . The next existence/well-posedness is usually called Lyons' universal 
limit theorem [Ij. 

Theorem 39 (Existence/well-posedness). 

a) Suppose the vector fields Vi are ('j—l)-Lipschitz and the vector field V is p-Lipschitz, 
for some < p < 1. Then the rough differential equation (jG.ip has a solution path. 

b) This solution path is unique if the Vi are ^-Lipschitz and V is (1 + p)-Lipschitz; it 
is in that case a continuous function of^. 
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We refer in the proof to propositions concerning Holder weak geometrc p- rough paths; 
using their analogue for branched rough paths proves the result in that case. Set a = 
min fl + -, -). 

Proof - a) Existence. Given e and t G [0,r], with ke s^ t < {k + l)e, set 

4 = iptke° Pke{k-l)e° ■ ■ ■ ° Peo){zo)- 

We prove that the paths (z|)o^t^T have a uniformly converging subsequence by showing 
that they form an equicontinuous family. 

First, for ke ^ s < t ^ {k + l)e, it is elementary to see that z| — z| = Ocjl^ — ■^k |, 
for some constant c independent of e. Now, given ^ s < t ^ T, with ke ^ s < 
{k + l)e, ie^t < {i + l)e, with A; + 1 ^ £, one has 



4 = Ptee° [ O Pin+l)ene] ° P(fc+l)e s(^s) • 



n=fc+l 

As every partition of any interval into sub-intervals of equal length is of special type 
I , one has by proposition [22] 

i-\ 

O P(n+l)ene = Ple{k+\y + Ol(|(A; - £)e\') 
n=A:+l 

provided t — s ^ 5. Since the maps pab are uniformly Lipschitz continuous, proposition 
[2T] gives 

zt = pts{zt) + OL+c{\t-s\''). 

The equicontinuity of the family of paths (2:^)o<e^i follows from the above inequality 

and the fact that Lots — Id ^ c[t — s[p. We see that any limit path satisfies (j6.2p 
for all ^ s ^ t ^ T with t — s ^ 6 hy sending e to 0. 

b) Uniqueness. To prove uniqueness of the solution when the Vi are 7-Lipschitz and 
V is (l-|-/9)-Lipscliitz, denote by z, any solution and remark that since it satisfies (16. 2p 
we have for any k ^ 1 

Zke - 4e = /'A:e(fc-l)£(^(fc-l)e) + Oc(e") - /3fce (fc-l)£ (^^(fc.i) J 

= Pke{k-l)eyP{k-l)e(k-2)e(Z(k-2)e) + Oc{e"-) j + Oc(e") " Pfce {fc-l)e(^(fc-l)e) 
= {Pke{k~l)€ ° P{k-l)e{k-2)e){Z(^k-2)€) + Oc(e") + Oc(e'') " Pk€ {k~l)e{zlk-l)e) 

as the pke(k-i)€ are Lipschitz continuous, uniformly with respect to k and e ^ 1. As all 
the pke(k-i)°' ■ '°P£e(l-i)e are Lipschitz continuous, with a uniform Lipschitz constant 
for any k ^ i and e ^ 1, by corollary [22] (for which we need the Vi to be 7-Lipschitz), 
a basic induction gives 

Zke-4e = {Pkeik^l) ° ■ ■ ■ ° Peo){zo) + Oc(^-j Oc{e'') - {pke{k-l) ° ' ' ' ° Peo) (zq) 

We conclude by sending e to along a subsequence for which z^ converges uniformly. 
> 
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Remarks 40. (1) The above notion of solution has a natural localized version. Given 

an open set U, a U-valued path (^;s)o^s<C with zq (^ U fixed and (^ ^ T, is said 
to solve the rough differential equation (|6.ip if there exists two positive constants 6 
and c depending on U such that (j6.2p holds for all ^ s ^ t < (^ with t — s ^ 6. 
The path z is called a local solution of (|6.ip . A solution {zs)os^s<( is said to 
he maximal if there is no other solution (^s)o<s<7 with C > C coinciding with 
z on [0,C]- The rough differential equation (j6.ip is said to have a solution if it 
has a maximal local solution for all bounded open sets. Replacing the boundedness 
hypothesis on V and the Vi by local boundedness, (|6.ip has a maximal solution path 
defined on some time interval [0,0 if the vector fields Vi are (7 — 1)-Lipschitz, but 
potentially unbounded; this solution path is unique if the Vi are j-Lipschitz, but 
potentially unbounded. By points (4) and (5) in remark [2S[ there is no explosion 
(C, = T) if the driving vector fields V, Vi have at most linear growth. See |11| and 
|19| for different proofs of that fact under slightly different hypotheses. 

(2) As in the preceeding section, the results of theorem [23 hold in an infinite dimen- 
sional setting, provided good tensor norms are used. It also holds in the Holder 
weak geometric case for the full rough differential equation. 

(3) Time reversal. Suppose the vector fields Vi are (potentially unbounded) j-Lipschitz. 
We see from remark\2^ that the time reversed solution path {zT-t)o^tsiT is the so- 
lution of the rough differential equation dz' = —Vdt-\-fX.. This simple remark can 
be used to construct martingales on a manifold whose final value is fixed. 

As we have only been working so far with Holder weak geometric rough paths 
rather than with rough paths associated with a general control uj{s,t) as can be 
done, we specialize this remark to that setting; there is no difficulty in extending it 
to the general framework. 

Let M be a d-dimensional manifold equipped with a connection V. Let J-iM.) 
stand for the frame bundle over M and Hi be the canonical horizontal vector fields 
on J^(M) associated with the connection V. Write vr : J^(M) — > M for the natural 
projection and H for {Hi, . . . , Hd). 

Definition 41. Let 2 < p < 3. An M-valued random process {xt)o<^t^T is said to be 
a special V-martingale if it is the projection of an T{M)-valued path {et)o^t!^T 
solving the rough differential equation de = HdX^, where X is the canonical Holder 

weak geometric p-rough path associated with an Mr-valued martingale X such that 

111- \ rt I I I- 

\Xt — Xs\ ^ c\t — s\p and J_, Xr ^ odXr\ ^ c\t — s\p . 

Given a final value xt in M and e^ G .^(M) with '7r(er) = xt, the time value 
of the solution in the above sense to the rough differential equation with unbounded 
vector fields de = H dji., if well-defined, provides a solution to the problem. If for 
instance M is isometrically C^ embedded in some W^ (which is not a restriction by 
Nash's theorem) the vector fields are the restriction to M of unbounded 'j-Lipschitz 
vector fields with at most linear growth, so the preceeding applies by point 1 above. 
Reparametrizing any W^-valued continuous (local) martingale by its p-variation up 
to the running time, it is possible to deal with the general case of a W -martingale 
in the setting of controlled rough paths described in point (4) of remark \15[ See 
|34| for a completely different approach based on backward stochastic differential 
equations. 
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7. Mean field stochastic rough differential equations 

We show in this section how the resuhs of sections [21 13] and [5] can be used to study some 
simple mean field stochastic rough differential equations. This kind of dynamics pops in 
naturally in the study of the large population limit of some classes of interacting random 
evolutions. The interaction holds through the dependence of the local characteristics of the 
random motion of each particle on the empirical measure of the whole family of particles. 
In a diffusion setting, each particle i would satisfy a stochastic differential equation of the 
form 

where fj,^ = i X]fc=i ^ W- ^ large industry has been devoted to showing that the limit 
distribution in paths space of a typical particle of the system when N tends to infinity has 
a dynamics of the form 

(7.1) dxt = b{xt, C{xt))dt + a{xt, C{xt))dBt, 

where C{xt) stands for the law of xt. Theorem |42] provides a well-posedness result for 
such a limit equation, in the context of rough differential equations. As emphasized in 
|35| . almost all the works in this area are set in the framework of a filtered probability 
space and rely crucially on some martingale arguments. On the other hand, the increasing 
importance of non-semi-martingale processes, like fractional Brownian motion, makes it 
desirable to have some more flexible tools to investigate equation (j7.ip in such contexts. 
The theory of rough paths developped above provides a nice framework for that. 

A few notations are needed to set the probem. Given 2<p<7^ [p] + l, we equip the 
set A1i(M ) of probability measures on M with the metric induced by its embedding in 
the dual of a (R"'): 

d(/i,z.) = sup {(g,/i) -{g,v)-ge C^W"), \\g\\^ ^ l}. 

This metric topology is stronger than the weak convergence topology. Given any posi- 
tive constant m, note that the set Lip(m-) of Lipschitz continuous paths from [0, T] to 
(A^ 1 (M'^) , d) , with Lipschitz constant no greater than m, is closed under the norm of 
uniform convergence, so it is a Banach space. 

Fix T > and suppose X is a random variable defined on some probability space 
(17, J-", P), with values in the set of Holder weak geometric p-rough paths over R^, on the 
time interval [0, T]; write X = 1©X^©- • -©XIp], and {J't)o-^t^T for the filtration generated 
by X. As above, F stands for a collection (Vi, . . . , Vi) of regular enough vector fields on 
M^. Given a Lipschitz continuous path "P, = (-Pi)o^ts£r in {Mi(^'^), d) and w G Jl, denote 
by x,{u)) the unique solution to the rough differential equation on paths 

(7.2) dxt = V{xt, Pt)dt + F{xt, Pt) dXt, 

where xq may be an integrable random variable independent of X. Denote by ^{V»)t 
the law of Xt] a solution x, to (17. 2p for which ^(V»)t = Pt is said to be a solution of 
the nonlinear rough differential equation (or mean field stochastic rough differential 
equation) 

dxt = V {xt, C{xt))dt + F [xt, C{xt)) dXt. 

Theorem[32]below provides conditions on V, F and X under which existence and uniqueness 
for solutions of this equation can be proved. We refer the reader to the forthcoming 
work |36] for a more general uniqueness statement. Note that the main assumption (j7.4p 
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below holds for instance for the Brownian rough path, the rough path above an Orsntein- 
Uhlenbeck process, or fractional Brownian motion with Hurst index no smaller than 2- 



Theorem 42. Let V,Vi,. . . ,Vi -.R"^ x Mi{R'^) -^ M'^ be vector fields on M^. Given any 
P G A^i(M ), we assume that V{-,P) is of class C^, with associated norm no greated than 
X, and that the vector fields Vi{-,P) are of class C'-^'^^, with associated norms uniformly 
bounded with respect to P, and that they satisfy the inequalities 



(7.3) 



max||Fi(-,P)-I^,(.,Q)|| v||F(-,P)-y(-,Q)|l ^Ad(P,g), 



i=l.I 



for all P,Q e Mi 



. Assume also that the random variables {X^'^ ) 



r=l..[p],/e[l,^I'- 



are 



integrable and that there exists for each ^ a ^ T, a positive random variable Ca, such 
that each of them satisfies the inequality 



(7.4) 



E 



X. 



r,I 
ba 



■J a 



^Ca{ 



for all ^ a ^b ^T, with sup ]E[Ca] < oo. Last, suppose that 



(7.5) 



Enixiri < 



oo. 



Then one can choose m big anough so that the map <l> has a fixed point in Lip(r?2). This 
fixed point is unique and depends continuously on the law of X if the vector fields Vi do 
not depend on their P-component. 

Proof - 1. Existence. We first prove that one can choose m big enough so that the 
map $ sends Lip(r7T,) to itself. Given such a path V,, the time-dependent bounded 
vector field V{x,Pt) is Lipschitz in its two arguments, so we can denote by ip^ the 
well-defined solution flow to the time non-homogeneous rough differential equation 
(rr2]l on R'^. We need to see that, for any function g G C'^{R'^), with ||5r||^ ^ 1, and 
any ^ s ^ t ^ T, we have 

{g,^r,)t) - {g,^r,),) ^ C{t - S), 

for some positive constant c ^ m, that is 



(7.6) 



E 



9{v>to(.^o)) - g{iPso{xo)) ^c{t 



s . 



Set Si = s + i2 "(t — s) and define the approximate solution flow 



(n),V 



^V 



V 



associated with the time non-homogeneous rough differential equation (17. 2p as in re- 
mark[25l(6): it converges uniformly to ip^ . The expectation E gi </?^(xo) ) —g( f'^oixo) 

is, by dominated convergence, the limit of E 






and writ- 



ing g\ pIq" [xq) ) — g\ PsQ (^o) ) as a telescopic sum, using the elementary identity 
(USD, gives 



E 



gi'^toi^o)] -aWsai^o) 



^v, 



A:=0 
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where we write z for 99^(xo), and p^ = p^^^k^i ° ' ' ' ° Psiso^ for 1 ^ A; ^ 2". With 
these notations, using the Taylor expansion of pf ^ given in (j3.7p . we see that 
^[gi^Toi^o)) - 9{^%{xo))] is the hmit of 



A;=0 



\P] 



r=i leliiY 

where o(-) does not depend on X. So, using (|7.4p and (j7.5p and the uniform bound- 
edness assumption in Cl^l"*"^ of the vector fields Vi{-,P), the quantity E[(7(99j^(xo)) — 
di^soi^o))] appears as no greater than the large n limit of 

2"-l [p] 



c(^-^)+EE E E[E[xr;;^3jj-,j(yK-,^sj5)(pr(^)) 

fc=0 r=l/G[l/|'- 



+ On(l) 



2"-l 



^ C(t - S) + C ^ E[C] (sfc+1 - Sk) + On{l) ^ C{t -S)+ On{l). 
fe=0 

Inequality (j7.6p follows by choosing m big enough. Fix Pq G A^i(R ) and define 
LipQ(m) as the subset of Lip(?Ti) paths with starting point Pq. As LipQ(?Ti-) is a convex 

compact subset of the set of continuous paths from[0,r] to (7Wi(M ),d), equipped 

with the norm of uniform convergence, Schauder's fixed point theorem applies and 
gives the existence of a fixed point of the map <&. 

2. Uniqueness. We suppose in this paragraph that the vector fields Vi{-,P) = 
Vi{-) do not depend on their A^i(M'^)-component. First, we prove that <l> is a strict 
contraction of LipQ(m), provided T is small enough. We use for that purpose the same 
telescopic decomposition as above, which is reminiscent of the well-known identity 
T( = Id + L hTrdr, satisfied by the semi-group {Tt)t^o of any Markov process with 
generator L. 

Fix a; € O and omit it in the notations of this paragraph. Denote by p^ the ap- 
proximate solution flow associated to the time non-homogeneous rough differential 
equation (17. 2p . with Qt instead of P(, for all ^ t ^ T. Given g G C^ijR. j, we see, 
using dominated convergence, that 

5(¥'?o(a^o)) -giv^Kii^o)) 



E 



lim E 

n— >oo 






is the limit of 



2"-l 

E 

A;=0 



E 



[{ 



9[P. 



Q 



.2 



J' 



rS2n^U + lS2"--k rS2n_f^S2n-k-l) 



9 [P. 



where p?„ 



P2"-k-l 



Q 

S2nS2n_i 
P 



Q 

' Ps2n-k + lS2"-k ° P- 



Q 

S2-^~k^2 



"-fc-1/ J 



-V 
° P2" 



-k-1 



(xo) 



S2n-k-lS2n-k-2 



■ o (Ol. o„, with the obvious convention concern- 
ing the summand for the first and last term of the sum. It follows from the time- 
inhomogeneous version of proposition [19] giving the Taylor expansion of Ps2Ti,-k+iS2-^-ki 
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that E 



5f(</?^(xo)) — g[(p^{xo)^ is no greater than the large n hmit of 



E 



2"-l 



V{;Q, 



'2"-k-iy\\oD 



2-"t 



+c||g|l^E[l+||Xp] 2^^". 



A:=0 

This upper bound is no greater than 

2"-l 

c \\g\\cf A 2-"t Y, d(^.„_,_i, Q.„_fe_0 + o„(l) ^ ct \\g\\c^ sup d(P„ Q,) + o„(l) 

since E[|[X|['''] is finite. As a result, we have for all t G [0, T] 

d(<I>(P.)t,<l>(Q.)t)= sup eU^^SMI -eUv^PoI^o)) 

llsllc^ssi 

^ cT maxd(Pt,(5t), 

iG[0,T] ^ ^ 

SO $ is a strict contraction provided cT < 1. 

As usual in the study of ordinary differential equations, the fact that <l> has a unique 
fixed point for any T follows from the fact that the above condition on T does not 
involve Pq = Qq. 



As the solution to the rough differential equation (|7.2p depends continuously on X, the 
map $, considered as a function of P, and of the parameter "law of X", is a continuous 
function of its two arguments. It is an elementary result that the unique fixed point 
of $ is then a continuous function of the parameter. > 



Remarks 43. (1) Following |35) . one can prove a propagation of chaos result for the 

particle system associated with (|7.2p . 

(2) Rough mean field games. Let us come back to the finite particle system de- 
scribed in heuristic terms in the introduction of this section. Suppose now that 
each individual can control his trajectory by tuning continuously a parameter u. 
Given a choice {v-t )o^t^T of control for the i^^ particle, its dynamics is described 
by the rough differential equation 



(7.7) 



dx 



(i) 



Si) 



,W 



V {x^ , /it , u\'> )dt + F {xf> ,fit,ur) dX 



(i) 



,{«)\ jy(*) 



Suppose now that each particle wishes to maximize some payoff of the form 



E 



J{s, Xs,fJ,s,Us) ds + C{xt) 



where the instantaneous cost J and the final cost C are given. Fixing the process 
(A'f )ot£s«;T and writing {tis)oKs^ 
tells us that the optimal payoff 



(/x^)o;gs^T and writing (/is)o^s!gT for it, the usual dynamic programming reasonning 



pt{x) ■.= supE / J{s,Xs,^^s,Us)ds + C{xT) 
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of a particle started at x at time t should heuristically satisfy the identity 
pt{x) = s\xp[J{t,x,^lt,u)e + ¥.[pt+e{xt+e)\ + o{e) 

= sup (j(t,x,^t,u)e + o(e) + dtpt{x) + E {Vpt){x, iJ.t,ut)e 

Yl Yl ^t+etiViPt)ix,Ht,ut) +Oc{ep 

r-=l/G[l,^r 

Under hypothesis (|7.4p , this equation amounts to a system of HJB equations with 
fixed terminal condition, whose form depends on how the expectations E[Xj'^^^] 
behave as e decreases to 0. If for instance 3 ^ p < 4 and the increments of X^ have 
zero mean, E[Xj_jf|,^] : 
we obtain the system 



zero mean, E[x'^^^^] = sf (e° + o(e")) for a < 1, andE[X^f^';] = ^j^'=(e + o(e)), 



dtpt{x) + sup( J{t,x,fj,t,u) + {Vpt){t,x,fj.t,u)+ Y. ^? (^j^j^fcf't)(i,a;,/xi,u)) = 0, 

u ^ i,j,k=l ' 



supl Y. Bl {VjVkPt){t,x,fit,u)) = 0. 



[1] 

[2] 
[3 

K 

[5] 
[6] 

[7 

[s; 

[9 



Suppose, in a given situation, that the HJB system is well-posed, with a maximum 
achieved at a unique point r{t,x, (/Us)t^s<T) for any continuous path (/^t)o^t<T- If 
each individual chooses this optimal control, equation (j7.7p specifies its dynamics. 
The mean field consistency condition couples the dynamics of a typical particle of 
the large N limit of the particle system with the above system of HJB equations by 
requiring that the law ut of Xf satisfies 

dxt = V\xt, ut,T{t, xt, (z^s)t^s^T) j dt + F\^xt, vt, r(t, xj, {vs)t^siiT)jdJLf 

Investigating existence and well-posedness for the coupled dynamics is a real chal- 
lenge; consult |37| for recent important advances on this area in the setting of 
Markovian dynamics driven by general Levy processes. 
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